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FINAL REPORT NCC1-213 (ODURF #152741) 


prepared by Hideaki Kaneko and Peter A. Padilla 


In the past twelve months (May 8, 1995 - May 8, 1996), under the cooperative agreement 
with Division of Multidisciplinary Optimization at NASA Langley, we have accomplished the 
following five projects: 

1. A note on the finite element method with singular basis functions 

2. Numerical quadrature for weakly singular integrals 

3. Superconvergence of degenerate kernel method 

4. Super convergence of the iterated collocation method for Hammersteion equations 

5. Singularity preserving Galerkin method for Hammerstein equations with logarithmic ker- 
nel. 

This final report consists of five papers describing these projects. Each project is preceeded 
by a brief abstract. The first author, Hideaki Kaneko, would like to thank Dr. Thomas Zang, 
Director of Division of Multidisciplinary Optimization at NASA Langley, for the financial support 
that he provided. 
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PROJECT 1 


HIDEAKI KANEKO and PETER A. PADILLA 

A Note on the Finite Element Method with Singular Basis Functions 

Recently there has been considerable interest in the finite element analysis that incorporates 
singular element functions. A need for introducing some singular elements as part of basis 
functions in certain finite element analysis arises out of the following considerations. The solution 
of certain problems, such as a field problem, exhibits highly singular behavior due to geometric 
features of the spatial domain. It is thought that an incorporation of singular elements that 
emulate the solution with the standard polynomial elements may perhaps be desirable. In order 
to make the computations of the finite element method with singular elements more efficient, 
Hughes and Akin ( The Finite Element Method, T.J.R. Hughes, Prentice Halt) established an 
algorithm for constructing interpolation functions that have the same interpolation properties 
of the Lagrange polynomials. 

We pointed out in this research that the aforementioned algorithm of Hughes and Akin is 
sensitive to the locations of the interpolation points that correspond to the singular basis func- 
tions. Specifically, we demonstrated numerically that the rate of convergence of a finite element 
solution varies according to the locations of these points. A general theoretical explanation is 
provided for this variance in the rates of convergence. 

A further invesigation is required toward the establishment of mathematical theory that 
guarantees the optimal rate of convergence of a finite element method that uses singular functions 
as part of its basis. 
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A Note on the Finite Element Method with 
Singular Basis Functions 

Hideaki Kaneko * 

Department of Mathematics and Statistics 
Old Dominion University 
Norfolk, Virginia 23529-0077 

Peter A. Padilla 

Multidisciplinary Design and Optimization Branch 
NASA-Langley Research Center 
MS-159 

Hampton, VA 23681 


Abstract 

In this note, we make a few comments concerning the paper of Hughes and Akin [3]. First 
it is noted that algorithm reported in [3] to produce a new set of functions is subject to the 
location of collocation points. Second we make a note concerning the rate of convergence of 
numerical solutions of the finite element method with singular basis functions. 
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1 Finite Element Method With Singular Basis Functions 


In the paper [3], Hughes and Akin made an interesting observation concerning the finite element 
analysis that incorporates singular element functions. A need for introducing some singular 
elements as part of basis functions in certain finite element analysis arises out of the following 
considerations. The solution of certain problems, such as a field problem [1], exhibits highly 
singular behavior due to geometric features of the spatial domain. On the other hand, in other 
circumstances, the solution is overwhelmingly affected by the nature of loading and the problem 
of singularity can be ignored. To satisfy both situations just described, it is thought that 
an incorporation of singular elements that emulate the solution with the standard polynomial 
elements may perhaps be desirable. This is the point that was exploited in [3] by Hughes and 
Akin . In order to make the computations of the finite element method with singular elements 
more efficient, they consider the following algorithm for constructing interpolation functions. 
A construction of such algorithm was motivated by the idea that “it is of practical interest to 
develop techniques for systematically defining shape functions for singularity modeling (and for 
developing special elements in general), which circumvent the interpolation problem” ([4] p.176). 
The algorithm that they developed go as follows: 


ALGORITHM Suppose that there are n shape functions N a , a = 1, 2, . . . , n which satisfy the 
interpolation property on the first m nodes r&, viz., N a (rb) = b a b, a,b = 1,2, . . ., m (m < n). 
Their idea here is to reshape N a * s so that the interpolation property is satisfied on all n nodes. 
The algorithm is given by 


Step 1 

Step 2 
Step 3 


N N m+ 1 (r)-ya=l^N m+1 (r a )N a (r) 

m+1 ■ N m+1 (r„)N Q (r m+1 ) 

N a (r) <r- N a (T) ^ N a (r m+1 )N m+ i{r), a - l,2,...,m 

If m + 1 < n, replace mbym + 1 and repeat Steps 1 to 3. 

If m 4- 1 = n, stop. 


To demonstrate this algorithm, we borrow one of the examples from [4], Let = 0, = 

and Tz = 1. The shape functions that we reconstruct are N\{r) — 1 - 2r, Ni{t) = 2 r and 
Ns(r) = r a where a representing some real number. Note that N a (rb ) = S a b, 1 < a, b < 2. An 
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application of the above algorithm gives 

Ni(r) *■ 
N 3 (r) * 
N 3 (r) * 


1 - 2r + I i-2(ij- 


2r — 2 


r°-2(i) a r 

1 - 2 ( 1 )° 


~ r n — 2(j’) a r 

. 1— 2(2 )° 


Of course, the newly defined shape functions satisfy 

N a (r b ) = 6 ab l<a,6<3. 


( 1 ) 


What is not addressed in [3] [4] is that the algorithm is subject to the location of the interpolation 
points rfc, m + 1 < b < n. Clearly, step 1 of algorithm does not work if these interpolation points 
are such that 

m 

N m +i(r m + 1 ) - N m +i(r a )N a (r m+1 ) = 0 ( 2 ) 

a=l 

Out of this observation, there seems to arise a profound and difficult problem in the area of 
approximation theory. The problem is important in that the success of the finite element method 
using the collocation scheme hinges on a resolution of this problem. To describe it, let W , * denote 
the Sobolev space, 

w} = {/|/W e T p (ft)} 

where /(*) denote the Arth generalized derivative and fi is a bounded region in R. The theory 
extends easily to higher dimensions. Now let U = span[i\T 0 ]"_ m+1 . Also denote an approximation 
space by S*. Here 5* is usually taken as the space of piecewise polynomials of degree k — 1 with 
length of each subinterval h. Our goal is to approximate each element of U © W£ by an element 
from 17 © S£ by interpolation. That is, if S£ = span[A a ]™_ 1 , then for each /€(/’© W£, we must 
find v £ U ® that satisfies 

v{n) = f{n) b = 1, ... .,n. (3) 

Denote the interpolation projector of U © W£ to U (& Sfc by Ph- Namely P/, is defined so that 

Phf(s) = v(s ) s € 0. (4) 

Notice that = Ph- In order to achieve a convergence by the collocation scheme in the finite 
element method, we must examine the following inequality. Here we assume that the number of 
singular basis functions, n — m, is fixed. 


( det [ JV a ( r b )] | ” 6=1 > e > 0 for all n. 


( 5 ) 
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This inequality is a necessary and sufficient condition for the algorithm of Hughes and Akin 
to work. It is important to remark at this point that the success of algorithm depends upon 
the existence of a solution to the interpolation problem (3) which in turn is equivalent to the 
condition (5). For each fixed index n (hence for fixed m ) , it is not difficult to find n — m 
interpolation points n,, m + 1 < b < n, that correspond to the singular basis elements N a , 
m + 1 < a < n for which the inequality in (5) is satisfied. What is difficult here is the question 
of locating n — m points for as many singular elements for which condition (5) is satisfied. The 
problem of finding n — m interpolation points for singular basis functions that work for all n 
is currently under investigations. In the following section, we proceed our discussion of the 
finite element analysis assuming condition (5). The analysis will provide information concerning 
the rate of convergence of numerical solutions. Interpolation examples at the end of the section 
show that rates of convergence are quite sensitive to locations of interpolation points for singular 
elements. 

2 Convergence Analysis 

When condition (5) is satistied, one can deduce the rate of convergence of the projector P ^ to 
the identity operator I. As is well known -e.g., [4], the convergence rate of such interpolation 
projectors determine the rate of convergence of the finite element method that uses collocation 
scheme. The following theorem of Cao and Xu [2] is useful. We sketch a proof for completeness. 

Lemma 2.1 Let X be a Banach space. Assume that U\ and U 2 are two subspaces of X with 
U\ C Vi- Moreover assume that Pi: X — * U\ and P 2 :X — + U 2 are linear operators. If P 2 is a 
projection, then 

Ik - P2x\\x < (1 + ||i 5 2||x)lk - Pix||x f° r al1 x e x. 

Proof: Let x 6 X. We write 

x - P 2 x = (x - Pix) + (Pix - P 2 x). 

Since P\x E Pi and U\ C U 2 , we have P 2 P\X = Pix. Hence, 

x - P 2 x = x - P\X + PvPiX - P 2 x 
= (I - P 2 )(x - P lX ). 


Ik - P 2 x\\x < (1 + ||P2||x)lk - P^x\\x for all X E X. 
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It follows that 



□ Using lemma and a well known fact about the approximation power of piecewise polynomials 
of degree k — 1 to approximate elements in Wf; [5], we see immediately that 

Theorem 2.2 Assume that y = u + v with u £ U and v £ W£. Then 

\\PhV - y\\ P < C7i n ||v (i) || 00 

where 1 < p < oo and C is a constant. 

To demonstrate the sensitivity of the location of interpolation points for singular elements, 
consider the following interpolation problem. 

EXAMPLE: Let f(x ) = yfx + %/l - x + x 2 . We wish to approximate / over [0, 1] by an 
element from UQS^, where U == spanf^/x, y/1 - xj. Let {x,}"_ 0 be the uniform partition of [0, 1] 
defined byzj = = and h = L, The interpolation points used to define an element 

from S'l are taken to be the zeros of the second degree Legendre polynomials transformed into 
[x;_i,x;] for i — 1,2, . . .,7i. The following data shows that (a) when the interpolation points for 
the singular elements are taken to be t\ — | and t^ —■ | for each n, the convergence is 0(h 1//2 ), 
whereas (b) when t\ = ^ and t? = 1 — ti, then the convergence is of the order 0(h 2 ). 


interpolation point t\ = 

1/2" 

1/5 

n = 4 

0.0513168 

0.0513168 

n = 6 


0.0235540 

convergence exponent = 

1.99 

0.56 


Table 1. Error and convergence rate data for the example 
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PROJECT 2 


HIDEAKI KANEKO and PETER A. PADILLA 

Numerical Quadratures for Weakly Singular Integrals 

In this report, a Fortran program for approximating weakly singular integrals is given. In 
particular, we are interested in approximating integrals of functions having algebraic or loga- 
rithmic end point singularities. Program is based upon the idea that was recently developed 
in the paper by Kaneko and Xu Gauss-type quadratures for weakly singular integrals and their 
application to Fredholm integral equations of the second kind- Mathematics of Computations, 
Vol. 62, (1994), 739-753. The idea of the quadrature scheme in the aforementioned paper is 
based upon the nonlinear splines approximation of weakly singular functions using a certain set 
of nonuniform knots that was originally developed by John Rice. 

The program developed here is useful for approximating solutions of weakly singular integral 
equations. This is demonstrated in Reports # 3 and # 4 in which the quadrature is used to 
approximate the solution of weakly singular Hammerstein equation that arises as a reformulation 
of the Dirichlet problem with a certain class of nonlinear boundary conditions. 
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Hideaki Kaneko’ 

Department of Mathematics and Statistics 
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Abstract 

In this report, Fortran programs for approximating weakly singular integrals are given. 
In particular, we are interested in approximating integrals of functions having algebraic or 
logarithmic end point singularities. Programs are based upon the idea that was recently 
developed in the paper by Kaneko and Xu [1]. 
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1 Theoretical Background 


In this report, we develop Fortran codes for Gauss- Legendre-type quadratures for weakly singular 
integrals. The idea of the quadratures was recently developed by Kaneko and Xu [1]. The integrals 
which we intend to approximate can be described in the form 

f f(x)w(x)dx (1) 

• Ja 

where / is a smooth function and w takes one of the following forms; 

1. (x — a)“(6 — x)P 

2. (x — a) a (b - x)Plog(x — a) 

3. (x — o)“(6 - x)^log(b - x) 

4. (x — a)“(6 — x)^log(x — a)log(b — x) 

where a,/3 > -1. To familiarize the reader with the basic theory that supports the current 
program, we shall consider the following problem. The reader who is interested in a more rigorous 
and broader development of the theory can find necessary materials in [1]. 

Consider 

1(f) = I f(x)w(x)dx (2) 

Jo 

where w(x) = x“, a > — 1, -i.e., w belongs to type 1 above and / is a smooth function. Define 

the index of parameter that measures the degree of singularity by 

2A; + 1 

9 ~ a + 1 

where k is a positive integer that corresponds to the degree k — 1 of polynomials used to approx- 
imate fw over [0, 1]. In the case of the logarithmic singularity w(x) = logx, we set a = 0. A 
partition 7r a of [0, 1] is defined according to the parameter a by 

n a : to = 0,h = n~ 9 ,tj =j v h, j = 2,3,...,n. (3) 

The more severe the singularity of w is at 0, the closer t\ is to 0. In the ordinary Gauss- Legendre 
quadrature scheme to approximate the integral JjJ f(x)dx , / is approximated by the polynomial 
Sk of degree k - 1 that interpolates / at the k points {x t '}^ =1 where {x,-}^ are k zeros of the 
kth degree Legendre polynomial transformed into [0,1]. The kth degree Legendre polynomial is 
described by 

2 X 
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The quadrature then takes the form 

jt, w if(xi) 

i=l 

where wf s are weights for the quadratue. It is well known [2] that the error of the Gauss- Legendre 
quadrature is given by 

/' H(. - if / € C< 2 ">[0, 1] (4) 

(2n)! Jo 

As evidenced by this formula for the error of the Gauss-Legendre quadrature, the scheme does 
not provide an optimal degree of approximation when the integrand / is not sufficiently smooth. 
It is exactly this reason why we require some special treatment for approximating the integrals 
in (2) that produces an optimal convergence rate. To achieve this goal, in the current method, a 
piecewise polynomial Sk of degree A: — 1 is constructed by the following rule; Sk(x) = 0,i6 [fo,ti) 
and Sk(x) is the Lagrange polynomial of degree k — 1 interpolating fw at for x 6 

i = 1,2, ... .,n-2 and for x £ [t n _i,t n ]. Here denote the zeros of the fcth degree Legendre 

polynomial transformed into Over the first interval, because of the endpoint singularity 

of fw , formula (4) is not at our disposal. Dispite this, we are able to recover the optimal rate of 
convergence due to sufficiently small size of the first interval that was determined by the index q 
of singularity corresponding to fw. To describe, in a more detailed way, the error of the current 
quadrature scheme, let 

rti+i rU + 1 

E k ,i= f(x)w{x)dx — / Sk(x)dx, i = 0, 1, ..., n - 1. (5) 

Jn Ju 

Then 

|£jfc,o| < ft„ | f{x)w(x) - Sk{x)\dx 

< /to \f{x)w(x)\dx 
— C Iff x a dx 

< Ct“ +1 = Cn ~ 2k ~ 1 . 

C denotes a generic constant whose value may change as it appears. For i > 1, by using (4), we 
obtain 

i^i < | (/ ~gii Wi>| i £+■ n,ii(x - 4>)^i 

< m a ~ 2k (t i+1 - u) 2k+i 
<ctr 2k (u+i-ti) 2k+i 

= C(i 9 t 1 ) Q - 2fc [(z + l ) 9 - z 9 ] 2 fc+ 1 tf +1 

< Cq 2k+1 (i + l)(9-l)(2fc+l)j9(a-2fc) f a+l 

= Cn~ 2k ~ 1 . 
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The last equality was obtained using 


q(a - 2k) = q(a + 1) - q{2k + 1) = (2 k + 1) - q{2k + 1) = (1 - q){2k + 1). 


Combining local errors above, the total error of approximation Ek is given by 

\E k \ = \J2E k ,i\<Cn- 2k ( 6 ) 

t=0 

obtaining the optimal convergence rate of approximation. 
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/* Author: Peter A. Padilla 


Summary: This program test the function sintegrate ability to integrate 

a function with weak singularities at both ends of the interval. 
The program calculates the integral in two ways, first we 
calculate the value by using sintegrate ability to automatically 
split the interval into two sections, in the second method we 
split the interval manually and feed sintegrate the appropriate 
parameters for each half of the interval. 


# include <stdio.h> 

# include <stdlib.h> 

# include <math.h> 

#include " sintegrate .h" 


/* Function Definitions 


double one (double x) 

{ 

return (1.0) ; 

} 

double fun (double x) 

{ 

return ( log (1-x) *sqrt(x) ) ; 

} 


/* Ma i n */ 

main(int argc, char argv[] ) 

{ 

double vl, v2 , v3 , v, (*gl)(), (*g2)(); 

int n; 

/* Example 3.2 of (1) */ 
n = 16; 
g2 = &one; 
gl = &fun; 

v = sintegrate (gl , g2 , 0.0, 1.0, n, 1, 1, 1, 0.5, 0.0); 
printf ( "Results : %g for n = %i\n",v,n); 

n = 24; 
g2 = &one; 
gl = &fun; 

v = sintegrate (gl , g2 , 0.0, 1.0, n, 1, 1, 1, 0.5, 0.0); 
printf ( "Results : %g for n = %i\n”,v,n); 

n = 28; 
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g2 = &one ; 


gl = &fun; 

v = s integrate (gl , g2 , 0.0, 1.0, n, 1, 1, 1, 0.5, 0.0) 
printf ( "Results: %g for n = %i\n”,v,n); 

exit ( EXIT_SUCCESS ) ; 



/ * Arguments : 


1) f is declared a pointer to a function that returns a double and takes 
a double as input, (double = double precision) . This pointer must 

be set to point to the weakly singular function for the interval in 
question, see test_driver. c for an example. 

2) weight is declared a pointer to a function that returns a double 
and takes a double as input. This pointer must be set to point to 
the smooth part of the kernel, see test_driver . c and test_driver_l .c 
for examples . 

3) a and b are the lower and upper bounds of the integration interval. 

4) n integer. Need I say more. 

’5) left_singularity is a flag, 0 if the singular point for the' interval 
is not at the left (point a) of the interval, 1 if it is. 

6) right_singularity is another flag. Can you guess how is defined? 

7) with_log is another flag. Set it to 1 if the function f is of types 
2, 3, or 4, i.e., it involves logO's. 

8) alpha is a double precision variable and must be set to the index of 
the singularity at x = a. 

9) beta is a double precision variable and must be set to the index of 
the singularity at x = b. 


*/ 


double sintegrate (double (*f ) (double x) , double (*weight) (double x) , 
double a, double b, int n, int left_singularity, 
int right_singularity, int with_log, double alpha, double beta) ; 
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/* Author: Peter A. Padilla 
Date: 9/21/95 

Summary: Implements the numerical quadrature algorithm for weakly singular 
integrals as described in the references . 

Ref: (1) H. Kaneko and Y. Xu, "Gauss-type quadrature for weakly singular 

integrals and their application to Fredholm integral equations of 
the second kind". Mathematics of Computations, Vol . 62, (1994), 

739-753. 

Calling symtax: See sintegrate.h 
Notes : 

The function sintegrate is the external interface to the code . No other 
function needs to be made accessible to the outside world. If other • 
interpolation polynomials are desired, then, the function legendre should 
be modified to implement the appropriate interpolation functions and the 
preprocessor constant K (line: "#define K 2.0") should be set- accordingly. 


*/ 

♦include <math.h> 
♦include <stdio.h> 
♦include <stdlib.h> 


/* Do not changed below this line unless you know what your doing 

♦define K 2.0 /* DO NOT CHANGE, K is determined by the 

interpolation function as described in the 
paper */ 


/* Function to generate a Partition */ 

double * points (int n, double a, double b, double q, int flag) 
{ int i; 
double tl; 
double *res; 

res = (double *) malloc (sizeof (double) * (n+1) ) ; 

if (q != 0) 
if (flag == 0) 

{ tl = (b-a)/pow(n, q) ; 

for ( i=0 ; i<=n; i++) res[i] = a + tl*pow(i, q) ; 

} 

else 

{ tl = (b-a) /pow(n, q) ,- 

for ( i = 0 ; i<=n,- i++) res[n-i] = b - tl*pow(i,q) ; 

} 

else 

{ printf ( "Error : q = 0 not allowed. \n" ) ; 
free (res); 
exit (0) ,- 

} 

return (res) ; 


/* Function to generate the interpolation points. */ 

double * inter_points (double p[] , int n) 

{ 

int i , j ; 
double eps ; 
double * tp; 



eps = l/sqrt(3.0); 

tp = (double *) malloc (sizeof (double) * (2*n) ) ; 

for (i=0 ; i<n; i++) 

{ j = 2*i; 

tp[j] = p[i] + (p [i+1] - p[i])/2.0*(1.0-eps); 
tp[j+l] = p[i] + (p[i+l] - p[i] )/2.0*(l. 0+eps) ; 

} 

return ( tp ) ; 

} 

/* Interpolation Function, Corresponds to k = 2 in the paper */ 

double legendre( double x, int j, int k, double * p) /* j=l..n */ 

{ • 
double point; 
double ret; 

point = (2 . 0*x - p[j] - p [ j -1 3 ) / (p [ j ] - p[j-l] ) ; 
if ((point < -1.0) || (point >= 1.0)) return (0.0); 
if (k == 1) ret = (1.0 - point*sqrt (3 ) ) /2 . 0; 
else if (k == 2) ret = (1.0 + point*sqrt (3 ) ) /2 . 0 ; 
else ret =0.0; 
return (ret) ; 

} 


/* Gaussian Quadrature with three points. */ 

double quadtdouble (*f)(), double a, double b, int j, int k, double * p) 

{ 

double cl , c2 , c3 , xl., x2 ,x3 , alpha, beta, Konst, temp, tempi; 
alpha = (b-a)/2.0; 
beta = (b+a)/2.0; 

Konst = alpha; 
xl = -sqrt ( 1 5.0) /5 . 0 ; 
x2 = 0.0; 
x3 = -xl ; 

cl = 1.0/ (3.0 *pow ( x3 , 2 ) ) ; 

c2 = 2. 0/3.0 * (3.0 - 1 . 0/pow(x3 , 2 ) ) ; 

c3 = cl; 

temp = Konst * (cl * (*f ) (alpha*xl+beta) *legendre (alpha*xl+beta, j , k,p) + 
c2 ' * (*f ) ( alpha *x2+beta) *legendre (alpha*x2+beta, j , k,p) + 
c3 * (*f)(alpha*x3+beta)*legendre(alpha*x3+beta,j,k,p)); 

return (temp) ; 

} 

/* sintegrate function */ 

double sintegrate (double (*f)(), double (‘weight) (), 
double a, double b, int n, int left_singularity , 
int right_singularity, int with_log, double alpha, double beta) 

{ 

double q_alpha, q_beta, .* grid, * gridl; 

double * interpolation_grid, * interpolation_gridl , * coefs, * coefsl; 
double * vals, value, * valsl, valuel; 
int i, j,k, index; 

if (with_log == 0) 

{ 

q_alpha = (2 . 0*K+1 . 0 ) / (alpha+1 . 0 ) ; 
q_beta = (2 . 0*K+1.0) / (beta+1 . 0) ; 

} 

else 

( 



if (alpha >= 0) 

q_alpha = 2.0*K+1.0; 

else 

q_alpha = (2 . 0*K+1 . 0 ) / (alpha+1 . 0 ) ; 

if (beta >= 0) 

q_beta = 2.0*K+1.0; 

else 

a beta = (2 . 0*K+1 . 0 ) / (beta+1. 0) ; 


} 


if ( lef t_singularity == 0) q_alpha =0.0; . 
if (right_singularity == 0) q_beta = 0.0; 

if ( (q_alpha != 0.0) && (q_beta == 0.0)) 

{ 

grid = points (n, a, b, q_alpha, right_singularity) ; 
in'terpolation_grid = inter_points (grid,n) ; 
coefs = (double *) malloc (sizeof (double) * ( 2*n) ) ; 
vals = (double *) malloc (sizeof (double) * (2*n) ) ; 

for { i=0 ; i<=2*n-l; i++) vals[i] = (* weight )( interpolation_gr id [i] ) ; 
for ( i = 1 ; i<=n; i++) 
for (k=l; k<=2; k++) 

{ 

index =2* (i-1) + k - 1; 
coefs [index] = 0.0; 
for ( j =0 ; j<=n-l; j++) 

coefs [index] = coefs [index] + quad(f, grid[j], grid[j+l] # i ( k, grid) 

} 

value =0.0; 

for (i=0; i<=2*n-l; i++) value = value + coefs[i]*vals[i]; 

free ((void *) coefs); 

free((void *) vals); 

free ((void *) grid); 

free ((void *) interpolation_grid) ; 

} 

else if ( (q_alpha == 0.0) && (q_beta != 0.0)) 

{ 

grid = points (n, a, b, q_beta, right_singularity) ; 
interpolation_grid = inter_points (grid,n) ; 
coefs = (double *) malloc (sizeof (double) * (2*n) ) ; 
vals = (double *) malloc (sizeof (double) * (2*n) ) ; 

for ( i = 0 ; i<=2*n-l ; i++) vals[i] = ( *weight) (interpolation_grid[i] ) ; 
for (i=l; i<=n; i++) 
for (k=l; k<=2; k++) 

{ 

index =2* (i - 1) + k - 1; 
coefs [index] = 0.0; 
for ( j =0 ; j<=n-l; j++) 

coefs [index] = coefs [index] + quad(f, grid[j], grid[ j+1] , i, k, grid) 

} 

value =0.0; 

for (i=0; i<=2*n-l; i++) value = value, + coef s [i] *vals [i] ; 

free((void *) coefs); 

free((void *) vals) ; 

free((void *) grid); 

free ((void *) interpolation_grid) ; 

} 

else if ( ( q_alpha ! = 0.0) && (q_beta != 0.0)) 

{ 

double mid; 
mid = (a+b) /2 . 0 ; 

grid = points (n, a, mid, q_alpha , 0 ) ; 
gridl = points (n, mid, b, q_beta, 1) ; 



interpolation_grid = inter_points (grid, n) ; 
interpolation_gridl = inter^points (gridl, n) ; 
coefs = (double *) malloc (sizeof (double) * (2*n) ) ; 
coefsl = (double *) malloc (sizeof (double) * (2*n) ) ; 
vals = (double *) malloc (sizeof (double) * (2*n) ) ; 
valsl = (double *) malloc (sizeof (double) * (2*n) ) ; 
for (i=0; i<=2*n-l; i++) 

{ vals[i] = (*weight) (interpolation_grid[i] ) ; 
valsl [i] = (*weiqht) (interpolation q ridlfi) ) 

} 

for (i=l; i<=n; i++) 
for (k=l; k<=2 ; k++) 

{ 

index =2*(i-l)+k-l; 
coefs [index] = 0.0; 
coefsl [index] = 0.0; 
for ( j =0 ; j<=n-l; j++) 

{ 

coefs[index] = coefs [index] +quad( f , grid[j ], grid[j+l] , i, k, grid) ; 
coefsl [index] = coefsl [index] +quad(f , gridl [ j ] , gridl [ j+1] , i, k, gridl) 

} 

) 

value = 0.0; 

valuel =0.0; 

for (i=0; i<=2*n-l; i++) 

{ value = value + coefs [i] *vals [i] ; 
valuel = valuel + coefsl [i] *valsl [i] ; 

} 

value += valuel; 

free((void *) coefs); 

free ((void *) vals); 

free ((void *) grid); 

free ((void *) interpolation„grid) ; 

free((void *) coefsl); 

free((void *) valsl); 

free((void *) gridl); 

free ((void *) interpolatio n g ridl ) ; 

) 

else if ( (q_alpha ==0.0) && (q_beta == 0.0)) 

{ 

printf ( "Internal Error: q == 0.\n"); 
exit(0) ; 

) 


return (value) ; 

} 



PROJECT 3 


HIDEAKI KANEKO and PETER A. PADILLA 

Superconvergence of Degenerate Kernel Method 

The degenerate kernel method is a classical method for finding approximate solutions of the 
second kind Fredholm integral equations ( x — Kx — f in operator form). The basic principle 
of the method is to approximate a kernel, a bivariate function, as a finite sum of univariate 
functions. An advantage of the method lies in its simplicity, whereas its disadvantage lies in the 
high cost of computations. If x n denotes a degenerate kernel approximation, then by the iterate 
of x n , we mean x' n = / + Kx n . 

We discovered that the rate of convergence of the iterates of degenerate kernel approxi- 
mations is determined by the method under which the kernel is decomposed. We proved and 
demonstrated numerically that, when the decomposition is done as a least squares approxima- 
tion or as an interpolation approximation using a certain set of interpolation points, then the 
iterates converge twice as fast as the original degenerate kernel solutions provided that the kernel 
is sufficiently smooth. 

To reduce the high computational cost, we propose to introduce a class of wavelets for a 
decomposition process. We expect to obtain a system of linear equations whose corresponding 
matrix is sparse as opposed to the normal full matrix that we encounter in the degenerate kernel 
method. 
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SUPERCONVERGENCE OF DEGENEREATE KERNEL MWETHOD 


H. Kaneko and P. Padilla 

Abstract 

In this paper, a general theory for the iterated operator approximation is developed. Some 
of the known results on superconvergence of various iterated schemes can be formulated as 
special cases of our theory. The method is'then subsequently used to obtain some results on 
the superconvergence of the iterated degenerate kernel method for the Fredholm equations 
of the second kind and the Hammerstein equations. 


1 Introduction 


(i.i) 


In this paper, we consider the Fredholm integral equation of the second kind 

x(t) — f k(t,s)x{s)ds = f(t) a<t<b 
Ja 

and the Hammerstein equation 

x(t) — f k(t, s)ip(s, x(s))ds = /(f) a <t < b 

Ja 

where k, f in (1.1) and (1.2) and ip in (1.2) are known functions and in both cases x is the function 
to be determined. A number of numerical methods are available for obtaining approximate 
solutions of (1.1) and (1.2). It is well known that popularly used methods such as the Galerkin 
method and the collocation method can be viewed as special cases of the projection method. 
The projection method for the Fredholm integral equation of the second kind (1.1) can be 
described as follows; Let X denote the Banach space of functions and {X n } a sequence of finite 
dimensional subspaces of X that is essentially dense in X,-i.e. cl U X n = X where clS denotes 
the closure of a set S. We denote a family of projections of X onto X n by {P n }. Now we define 

f b 

Kx(t)= I k(t,s)x(s)ds for each t € (a,b). (1.3) 

Ja 

Equation (1.1) can be written in operator form as 


x — f + Kx. (1-4) 

Here x is assumed to be an element of X. The underlying principle that supports the projection 
method is to seek an element x n G X n for which the residual r n = / - (x n — Kx n ) disappears 
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under the projection P n , -i.e. P n r n = 0 for each n. In the Galerkin method, the projection is 
orthogonal. On the other hand, in the collocation method, the projection is interpolatory. A 
similar discussion for the projection method for the nonlinear Hammerstein equation (1.2) will 
be left to the reader (see, e.g. [2] [3] [4] [8] [9]). 

Let x n denote a numerical solution of equation (1.4). The (Sloan) iterated solution x' n 
corresponding to x n is defined by 

*n(0 = /(*)+/ k(t,3)x n (s)ds. (1.5) 

Ja 

In the Galerkin and the collocation solutions x n , it is known [6] that, under suitable conditions 
on the kernel k and the forcing term /, their corresponding iterated solutions converge to x more 
rapidly than x n does to x, a phenomenon commonly known as superconvergence. One of the 
purposes of this paper is to develop a general iterated approximation scheme for equation (1.1). 
This will be done in Section 2. The generality of the theory is demonstrated by showing that 
the results on the superconvergence of the Galerkin and the collocation iterates examined by 
Sloan et al are special cases of our iterated operator approximation scheme. We note that the 
results obtained in [6] concerning the superconvergence of the iterated Galerkin and the iterated 
collocation method have been recently generalized to hold for Hammerstein equation (1.2) [10] 
[7]. Moreover, in these papers superconvergence of the iterated solutions for weakly singular 
Hammerstein equations are also observed. It appears that these results on superconvergence of 
the weakly singular equations are new even for the Fredholm equations. 

Results obtained in Section 2 are subsequently used to establish the superconvergence of the 
iterated degenerate kernel method for equations (1.1) and (1.2). The degenerate kernel method 
for solving equation (1.1) is discussed in [1] and that for equation (1.2) is discussed in [11]. It 
turns out that the superconvergence of the iterated degenerate kernel method depends’ upon 
the ways in which the kernel k in equations (1.1) and (1.2) are decomposed as a finite sum of 
products of univariate functions. This will be illustrated in detail in Section 3. In Section 4, the 
iterated degenerate kernel method for Hammerstein equation (1.2) is discussed. Even though the 
material in Section 4 does not have direct relation with the general theory presented in Section 
1, we feel that it is appropriate to include it here since it extends the results of Section 3 as well 
as those in [11]. Also included in Section 4 is a discussion on a computational consideration for 
solving the system of nonlinear equations that must be solved in order to obtain the degenerate 
kernel solution of (1.2). We make use of the devise introduced by Kumar and Sloan [12] to 
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‘linearlize’ the system of nonlinear equations. That is, the integrals that must be evaluated 
repeatedly at each stage of the iterations, when a conventional iterated scheme is employed, are 
now evaluated only once at the beginning of computation. 

Examples are provided in Sections 3 and 4 to demonstrate our theory. 


2 The Iterated Approximation for the Fredholm Integral Equa- 
tions of the Second Kind 

In this section, we consider the Fredholm integral equation of the second kind given by (1.1). 
We denote by C[a, 6] the Banach space of all continuous functions defined on [a, 6] equipped 
with the uniform norm || •Hoc. Also we denote by L p [a, 6], 1 < p < oo, the Banach space of 
measurable functions whose pth power are integrable (1 < p < oo) or the space of essentially 
bounded functions ( p — oo) on [a, 6]. In equation (1.1), we assume that 

feC[a,b\ (2.1) 

and that, with k t (s ) = k(t,s), 

f b 

lim / \k t (s) — k T (s)\ds = 0 for each r <= [a, b]. (2.2) 

t~*T Ja 

Under the condition (2.2), the linear operator K defined by (1.3) is a compact linear operator 
of C[a,b\ into C[a,b]. In order to establish a general iterated approximation scheme, we assume 
that {K n } is a sequence of operators converging to K in operator norm. That is, 

\\K n — K\\ p — ► 0 as n — * oo for some 1 < p < oo. (2-3) 

For each n > 1, we assume that we have an equation whose solution approximates the solution 
of (1.4). We denote this approximating equation by 

x n = fn + K n x n . (2.4) 

For example, in the case of the projection method, equation (2.4) is identified by letting K n = 
P n K and f n = P n f where P n is a projection of C[g, 6] onto a finite dimensional subspace of 
C[a, 6]. In the case of the degenerate kernel method, K n denotes the finite rank separable 
operator, -i.e. K n x{t) = fa J2?=i o-i(t)bi(s)x(s)ds where {o,}" =1 is a linearly independent family 
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of functions and /„ = / for each n > 1. We define the iterated approximation corresponding to 
(2-4) by 

x' n = f + Kx n . (2.5) 

As was indicated in Introduction, the iterated approximations for the Galerkin and for the 
collocation methods exhibit, under suitable smoothness conditions on the kernel k and on the 
forcing term /, the global superconvergence. It is not known, however, that the similar super- 
convergence result can be obtained for the iterated approximations for the degenerate kernel 
method. We are now in a position to prove the main theorem of this section. Advantages of the 
theory presented are twofold. First, it can be used to describe the superconvergence results of 
Sloan et al mentioned as special cases of this theorem. Second, it can be used to establish the 
superconvergence of the iterated degenerate kernel method. The latter will be done in Section 
3. 

Theorem 2.1 Consider equation (1.1) in a Banach space (X, || - ||). We assume that 1 is not 
an eigenvalue of K. Let x n and x' n satisfy equations (2.4) and (2.5) respectively. Then, for 
sufficiently large n, there exists a constant c > 0, independent of n, such that 

11* - <11 < c{\\K - K n || 2 + || K(K- K n )x || + || K - K n \\ || / - / n || + ||A'(/ - /„)||}. (2.6) 

Proof: From (1.4) and (2.5), 

x-x' n = K(x - x n ). (2.7) 

Applying K on both sides of (1.4) and (2.4), we obtain 

Kx = Kf + K 2 x (2.8) 

and 

Kx n = KK n x n + Kf n . (2.9) 

It follows from (2.8) and (2.9) that 

K(x-x n ) = K 2 x-KK n x n + K(f-f n ) 

= K(Kx - K n x n ) + K n (Kx - Kx n ) - K n (Kx - Kx n ) + I((f - f n ) 

= K n (Kx - Kx n ) + (K - K n )(Kx - Kx n ) + K(K - K n )x n + K(f - f n ). 

(2.10) 
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Since \\K n - K\\ -*Oasn-»oo and (I - K) 1 exists by assumption, we conclude [1] that 
(/ — A' n ) _1 exists and uniformly bounded for sufficiently large n. Therefore, 


K(x - x n ) = (I- KJ-'iiK - K n )(Kx - Kx n ) + K(K - K n )x n + K(f - /„)}. 


Taking the norm on both sides, 

\\K(x - x„)|| < ||(/ - A' n )- 1 ||{||A' - K n \\\\K\\\\x - * n || + \\K(K - K n )x n \\ 

+\\K(K - K n )(x - x n )\\ + || K(f - fn)\\} 

< ||(/ - K n )-'\\{2\\K\\\\K - tf n ||||* - x n || + || K(K - K n )x n 
+\\K(f-f n)||}- 


Since 


x x n — Ax A„x n -I - f f n 

— Ax K n X *1“ K n X ^n*^n T f fn 


we obtain 


(/ - A'„)(x - x n ) = Kx - K n x + /-/„, 


( 2 . 11 ) 


x - x„ = (7 - K n )~ l {(K - K n )x + f- /„}. (2.12) 

From (2.11) and (2.12), 

II*- <11 =ll#(*-*n)ll 

< c{||iT - K n ||||x - x n || + \\K(K - K n )x n \\ + \\K(f - f n ) ||} 

< c{\\K - K n \\ 2 + ||A - K n \\\\f - /n|| + \\K(K - K n )x n \\ 

This completes our proof. □ 

Corollary 2.2 For the iterated approximation scheme (2.5), if /„ = / for all n in (2.4), then 


II* ~ <11 < <11 K - K n f + || K(K - AT n )*n||}. 

In the next section, this corollary will be used to establish a superconvergence result for 
the iterated degenerate kernel scheme. In order to see that Theorem 2.1 includes the results of 
superconvegence of the iterated Galerkin and the iterated collocation schemes of Sloan et al [6], 
we need some definitions. First we let W™, 1 < p < oo, m nonnegative integer, denote the 
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Sobolev space of functions defined over [a, 6]. Namely g G W™ if and only if gW G L p [a, 6] for 
k = 0, 1 , . . m where g M denotes the A:th distributional derivative of g. Let 


II„ : a = t 0 < *i < • • - < t n = b 

be a partition of [a, b]. Let h = maxi<j< n (t,- — ti_j) and assume h 0 as n -* oo. Let u be 
an integer and r a positive integer such that 0 < v < r. Let S^ n denote the space of splines 
of order r and continuity is, namely <p n G S^ n if and only if ip n is a piecewise polynomial of 
degree < r — 1 on each X{] and has v — 1 continuous derivatives on (a, b). Let P^ denote 
an orthogonal projection of C[a,b] onto S^ n . In the Galerkin method, equation (2.4) becomes 

x°-P?KxZ = PZf ( 2 . 13 ) 


-i.e. K n = P®K and f n = P^ f. The corresponding (Sloan) iteration approximation to (2.5) is 
given by 

= / + **?. (2.14) 

If / G Wp,(m > 0), then there exists tjj n G 5£n (0 < v < r) such that 

(2.15) 


where c is a constant independent of n and ||/|| m>p = o 11/^ lip ( see e -g- [13]). Under the 

assumption of the quasiuniform mesh, -i.e. 


it can be shown that 


h 

mini<j<n(L' f»-i ) 


< c 


for each n, 


SU P \\Pn\\L p -*L p < C 2 - 

n 


for some constant c-i. Since 


(2.16) 


\\f-P?f\\r =\\f-^n + P^n-PSf\\p 

<(l+l|ff 11)11/ -^n||p, 

from (2.15), (2.16) and (2.17), we obtain 


(2.17) 


ll/-P n G /llp<^ min{m ’ r} ||/lk P . 


(2.18) 
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Now let £(t) = f* k(t,s)z%(s)ds. Then 




= I £*(«, «){{(«) - 

= 1(*« -Vn,(- p% f )l for every y>„ e 

<iiA.-v>~iu?-ffeiu , 


(2.19) 


where L + L = 1 with convention that if p = 1, then q = oo. In (2.19), we have used the 
orthogonality in the third equality and the Holder inequality in the last step. If k t G W™, then 
from (2.15) there exists <p n € S" n such that ||fc t - <p n \\ q < c / l imn { m - r }||fc[ m )|| m9 . Finally from 
(2.19) we obtain 


\\K{K - KJxZWoc < ch 2min { m > r }. 


Similarly, we can show that whenever / € W™, 


\\K(f - Ff /)||oo < ch 2uda{m ’ r} 


and that 

\\K - FTnlloo < Ch m . 

Using the estimate (2.6), the above discussion leads us to the following corollary. 

Corollary 2.3 (Graham, Joe and Sloan [ [6] : Theorem 4-1]) Let and x„' denote the solutions 

for (2.13) and (2.14) respectively. Suppose that k t € W™ (0 < m < r) with ||fci|| m , g bounded 
independently oft and that /,£ € W™ where £(t) = k(t, s)x^(s)ds with x„. Then 

where c is independent of n. 

Now in the partition II n , for each i, we select {Uj} T j = i such that 

Xi-l < til < til < • ■ ■ < tir < Xi. 

Let denote the interpolatory projector of C[a,6] onto S(. n defined by P£x(tij) = x(Uj) for 
each i = 1, . . . , n and j — 1,2 , . . . , r . In the collocation method, equation (2.4) becomes 

x n ~ Pn Kin = P% f (2.20) 
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-i.e. K n = P£ K and f n = P% f. The corresponding iterated collocation solution is defined by 


*? = / + ffa£. (2.21) 

As in Corollary 2.3 of the iterated Galerkin method, to see that the iterated collocation 
method of (2.21) is a special case of Theorem 2.1, we must examine the terms in the right side of 
(2.6). The second term of (2.6) in this case is analyzed as follows: Let y{t) = / a b k(t, s)i£(s)ds. 
Then 

K(K - K n )xC (t) = (k t , y - P%y) 

(kt - y - PSv) + On,*. (I - Pn){y ~ V’n)) ( 2 - 22 ) 

+(tpn,t,(I - P°)i> n ). 

where y n< t € S^ n and ip n e Sf n . Now arguing exactly as in the proof of theorem 4.2 [6], we 
obtain 

P(tf - K n) x n\\o° < c/l 2min ^’ m+r} 

where c is a constant independent of n. The other terms in (2.6) can be bounded similarly. 

Corollary 2.4 (Graham, Joe and Sloan [[6] : Theorem 4-%]) Let x n an< ^ x n be the solutions 
of (2.20) and (2.21) respectively. Suppose f € C[a,b], x € W{ (0 < l < 2 r) and kt 6 WJ 71 
(0 < m < r), with ||A:t|| m ,i bounded independently oft. Then 

||* - *?' lloo < 

where c is independent of n. 

3 The Iterated Degenerate Kernel Methods 

The purpose of this section is to use Theorem 2.1, Corollary 2.2 in particular, to establish 
superconvergence of the iterated degenerate kernel method. A generalization to the iterated 
degenerate kernel method for the Hammerstein equations will be done in Section 4. 

• Consider equation (1.1). The degenerate kernel method for approximating the solution of 
(1.1) requires us to approximate the kernel k in (1.1) by some degenerate kernel whose general 
form can be described as 

n n 

k n {s, t) = y ] ajj<pi(s)(pj(t) (3.1) 

i=i j= l 
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where {¥>,}" = i is a set of linearly independent functions in C[a,b], The operator K in (1.3) is 
then approximated by a sequence of operators 


K n y{t) = [ k n (t,s)y(s)ds. 

Ja 

Subsequently an approximate solution x n is found by solving 

rb 

x n (t) - / k n (t,s)x n (s)ds = f(t) a <t <b. 

J a 

Equation (3.3) can be written as 

*n(0 - EwW{E I aij<Pj(s)x n (s)ds} = /(/) a <t <b. 

■+ i Ja 

i=l j=l 

If we put 

rb 

Ci = J2 / aim{s)x n (s)ds, 
i= l Ja 


then x n can be written as 


®n(0 = + 


i = 1 


(3.2) 


(3.3) 


(3.4) 


(3.5) 


Upon substituting (3.5) into (3.4), we obtain the following n x n system of linear equations for 


Ci. 


n « rb n rb 

aij<pj(s)(p k (s)ds = Y^ a im( s )f( s ) ds 1 < i < n. (3.6) 

z=i j=i Ja j=i Ja - 

Finally, once these c,-’s axe found by solving (3.6), equation (3.5) gives the required approximate 
solution for the degenerate kernel method. Equation (3.3) is written in operator form as 


x n - K n x n = f (3.7) 

which is a special form of (2.4) with f n = f for all n. When the degenerate kernel solution 
x n is iterated as in (2.5), an interesting question arises. The question is of course under what 
conditions superconvergence of the iterates in this case is guaranteed. It is the purpose of this 
section that we provide some answers to this question. It turns out that the superconvergence 
of the degenerate kernel method hinges critically upon the ways that one decomposes the kernel 
k in (3.1). Here we demonstrate two different methods that guarantee the superconvergence of 
the iterates of the degenerate kernel method. 

In the first method, we examine the least-squares approximation. For each positive integer 
k , assume that a partition 11* satisfies the quasiuniform condition (ref. Section 2) 

— y < c for each k 

mini <j<n(^t l ) 
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Let n denote the dimension of the spline space Sf. k and Bi,B 2 ,...,B n be the H-spline basis 
for S" k . As in Section 2, r and v are integers such that 0 < v < r. Assume that k n (t,s) is the 
least-squares approximation of k(t,s),- i.e. assume that a,j in (3.1) are such that 

/ / !*(*, o.ijBi{s)Bj{t)\ 2 dsdt = min f i | k(t, s)-J2J2 hjB^Bj^dsdt. 

Ja Ja i=l j=l o.jenJa J a i=l j=l 

We are now in a position to present the first theorem of this section. 

Theorem 3.1 Let x be the solution of (1.1) and x n the solution of (3.7) where (3.1) is defined 
by the least-squares approximation. Assume that k £ W™([a, b] x [a, 6]), 0 < m < r, and 
k s (-)x(-) £ W 2 ([a,&] x [a, 6])) for each s £ [o, b], where 0 < l < r and k s (t) = k(s,t). Then 


II* - *»lloo = 0{h*) 


with v = min{m + l, 2m}. 


Proof: Using corollary 2.2, we obtain 

= 0(h? m ) + 0(\\K(K - if„)x n || 00 ). 


x — x 


nlloo 


(3.8) 


Hence we only need to estimate the order of convergence of \\K(K — ivnJinHoo. Note that 

I K(K - K n )x n (t ) | = | / a 6 k(t, u) J?lk(u, s) - k n (u, s)]x n (s)dscfa| 

= fa fa »)[*(«> s ) - fc n(«, s)]x n (s)dsdu\. 

Let i/j t (u,s) = k(t, u)x n (s) and let cp n (u,s) = X)?=i X l]=i bijBi(u)Bj(s) be any element that is a 
tensor product of H-splines. Then since k n is the least-squares approximation of k, 


n ip n (u, s)[k(u , s) — k n (u , s)]dsdu = 0, 

_ 


therefore 


\K(K - AT n )x n (f)| = | f f [V>t(u, s) - <p n (u, s)p(ti, s) - k n (u, s)]dsdu. 
Ja Ja 

Applying the Cauchy-Schwartz inequality, 

\K(K - K n )x(t) | < - <pnh \\k - k n \\ 2 . 


Here of course jj • ||2 denotes the L 2 norm defined on the space of bivariate functions W^fa, b] x 
[a, 6]). Noting that \\k — fc n || 2 = 0(h m ) and \\ip — ip n \\2 = 0(h l ), (3.8) proves the desired result. 
□ 
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The second method that produces superconvergence of the iterates of the degenerate kernel 
solutions is based upon approximating k by interpolation. Let 


II/v : a — to < t\ <•■•<. t/v = b 

be a partition of [a, 6]. Let h = maxi<;<jv(f; - L'-i ) and assume as in Section 2 that h -* 0 as 
N —> oo. Let £i,£ 2 , . . .,£ r be the zeros of the rth degree Legendre polynomial in [-1,1]. We 
shift these points to each subinterval i = 1,2, ...,N to obtain Denote the 

interpolation polynomials by (pij, -i.e. 


<Pij(T a ,p) = 

An approximating degenerate kernel k n (n = Nr) is now defined by 

N r N r 

k n(s, <) = IXEE *( T iji T ctj)vM<p a j3(t). ( 3 - 10 ) 

t=l jfssl a=l /3=1 


1 if (i,j) = {a, /3) 
0 if (i,j)^ (a, /3) 


(3.9) 


Let the interpolation projector of C([a,6] x [a, 6]) into the tensor product space 5° n ® S°„ be 
denoted by P n . That is, 


P n k(s,t) = k n (s,t) 


where k n is defined in (3.10). Also let x£ denote the solution of equation (2.4) when K n possesses 
the kernel k n defined by (3.10) and x% denote the corresponding iterate defined by (2.5). The 
following theorem demonstrates the superconvergence of the iterated degenerate kernel method 
when the kernel is decomposed by interpolation. 


Theorem 3.2 Assume that in equation (l.l)k G W 1 m ([a, 6]x[a,6]) ; 0 <m<r, and k s (-)x n (-) G 
VPj([a,6] x [a, 6]), 0 < / < 2r, where kt(u ) = k(t,u). Then 

||x - x^Hoo = 0(h J/ ), v = min{2m,r + /}. 

Proof: As in the proof of Theorem 3.1, we need to estimate the error of || K(K — -Kn)x n ||oo. 
By taking <p n G Sf n ® S^ n and ip n G S° ® 5° >n , for each t G [a, 6], 

K(K- K n )x(t) = / a 6 fc(f,u)/ a 6 [fc(n,s) - fc n (u, s)]x n (s)dsdu 
= Ja fa Wi u ) x n{s)[k(u, s) - k n (u , s))dsdu 
= (fc t («)x„(s), k(u, s ) - k n (u, s)) 

= (fc«(u)a:n(s) - <p n (u, s),k(u, s) - k n (u, s)) 

+(ip n (u, s), (I - P n )(k(u, s) - i) n {u, s)) + (< p n (u , s), (I - P n )^n(«, «)) 
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The rest of proof follows by an argument similar to the one give in the proof of theorem 4.2 of 
Graham, Joe and Sloan [6]. A modification needed here however is to accomodate the bivariate 
functions. An approximation of a bivariate function by an element from a space of tensor product 
of finite dimensional univariate functions and its approximation power is well ducumented, -e.g., 
see [5]. □ 

4 The Iterated Degenerate Kernel Method for Hammerstein 
Equations 

In this section, we extend the results of the previous section to obtain superconvergence of the 
iterated degenerate kernel method for Hammerstein equation described in (1.2). We assume 
throughout this section, unless stated otherwise, the following conditions on fc, / and ip: 

1. lim ( — jj/uj fcr||oo = 0, t G [fl, 6], 

2. M = sup a <,,< 6 / a 6 \k(t,s)\dt < oo; 

3. / € C[a, 6]; 

4. ip{t,x) is continuous in f € [0,1] and Lipschitz continuous in x G (- 00 , 00 ), i.e., there 
exists a constant C\ > 0 for which 

\ip(t,xi) - ip(t,x 2 )| < Cil^i - x 2 \, for all x x ,x 2 G (- 00 , 00 ); 

5. the partial derivative ip( a,b ) of ip with respect to the second variable exists and is Lipschitz 
continuous, i.e., there exists a constant C 2 > 0 such that 

| ip (a,b \t, Xi) - ^ (0,1) (t, x 2 )\ < C 2 \xi - x 2 \, for all ®i, x 2 G (- 00 , 00 ); (1.2) 

6. for x € C[0, 1], ip(., x(.)), ^ (0 ’ 1) (., x(.)) G C[a, b ]. 

Additional assumptions will be given later as needed. A comprehensive study of the degenerate 
kernel method for Hammerstein equations was made by Kaneko and Xu [11]. The degenerate 
kernel method for Hammerstein equations consists of replacing K in (1.2) by k n of (3.1) and 
approximating the solution x of equation (1.2) by x n which satisfies 

f b 

Xn (t) - / k n (t,s)ip(s,x n (s))ds = f(t ) a < t < b. (4.1) 
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Following analogously the development made in (3.5) and (3.6), for 

" rb 


x n can be written as 


74 rt> 

= a ij'W( a )lK a » *"(«))*, 

j=l Ja 


Xn (*) = /(O + ^CiWW- 
«=1 


(4.2) 


(4.3) 


Substituting (4.3) into (4.2), we obtain the following n nonlinear equations in n unknowns 
’ rCn? 


"i = 22 I f(s) + Zu c ‘W( s )) ds > 1 < i < n. 

j = l Ja l=i 

KVx(t) = f k(t, s)0(s, x(s))ds 
Ja 


Define 


(4.4) 


so that (1.2) can be written as 

x — K9x = f. 

Similarly we write equation (4.1) as 

x n ~ K n ^Sx n — f 

The iterated solution x' n is now obtained by 

x' n = f + K*x n . 


(4.5) 


(4.6) 


(4.7) 


The Frechet derivative of K9 at tp 0 is denoted and defined by 

(fir$)'(<po)(¥0(*) = / k(t,s)^ 2 (s,<po(s))(p(s)ds 
Ja 

with i / , 2 denoting the first partial derivative of 0 with respect to the second variable. The 
following theorem describes superconvergence phenomenon of x' n to x. Here we assume that, 
the decomposition (3.1) is done via interpolation (the second method described in the previous 
section). The case of the least-squares approximation is similar and will be left to the reader to 
supply the detail. 

Theorem 4.1 Assume that in equation (1.2) k € IF™ ([a, 6]x[a, 6]), 0 < m < r, and fc s (-)0(-, x„(-)) € 
VF(([a, 6] x [a, 6]), 0 < / < 2r, where x n is the solution of (4-6). Assume also that 1 is not an 
eigenvalue (if0)'(£) for each £ between f + Kn'Sxn and f + K'Sx. Then 

ll x “ x n lloo = 0(h"), V - min{2m,r + /}. 
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Proof: From (4.5) and (4.6), 

x - x' n = KVx - KVx n . (4.8) 

Now 

KVx - KVx n = K*U + KV x) - KV{f + K n 9x n ) 

= (KV)'(6(f + K n Vx n ) + (1 - 0)(f + KVx))(K*x - A' n $i„) 
for someO < 9 < 1 

= Kg{KVx - K n *x n + (KVx - KVx n ) - (K9z - K$x n )), 

where Kg = (Kty)'(8(f + K n ^x n ) + (1 — 0)(f + K9x)). Since 1 is not an eigenvalue of Kg for 
each 0 < 6 < 1, we obtain 

KVx - KVx n = (/ - Kg)~ 1 Kg(K - K n )Vx n . (4.9) 

Combining (4.8) and (4.9), and taking the norm on both sides, we obtain 

II* - *'Joo < c\\K e (K - K n )V I„||oo. 

Arguing as in the proof of Theorem 3.2, we obtain the desired result. □ 

Finally we consider a computational problem associated with (4.4). It is customary that 
the system of nonlinear equations (4.4) is solved by an iterative scheme. For example, the fixed 
point iteration scheme for (4.4) is to generate {c[ fc ^}" =1 for k > 1 by 

c .' fc+1) = £/ J2c\ k) (pi(s))ds, 1 <*■<«. (4.10) 

j=l Ja 1=1 

At each step of iteration, the integrals in (4.10) must be computed since the integrands contain 
the different values of c\ k \ To circumvent this difficulty, we propose the following device that 
is originated from [12]. We let 

z n (t) = i{)(t,x n (t)) (4.11) 

where x n is defined in (4.3). We have, assuming that k n takes the form of (3.1), 

z n {t) - 4>(t, f(t) + f ¥>j(s)z n (s)ds). (4.12) 

«=i Ja j= l 

The equation (4.12) can be solved by the collocation-type scheme that was developed by Kumar 
and Sloan [12]. Namely let {rji}? = i be n functions defined on [a, 6] and let {tj} r - =l be n distinct 
points for which 

detfafo)) ± 0. (4.13) 
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z n in (4.11) is now approximated in the form a jVj- The afs can be found by solving the 
following nonlinear equations. Note that the constants ctj’s are moved out of the integrals. This 
makes the repeated computations of the integrals unnecessary when the system of nonlinear 
equations is to be solved by an iterated scheme. 

]£ a jVj(tk) = ^(<fc, /(*t) + ]C a H Vi(*k) Y, ai [ it, V'MVii*) 11 *)- ( 4 T4) 

j = 1 i=l 1=1 Ja j=l 

for 1 < k < n. If we denote A = [qj(ti)j and the right side of (4.14) by ipi(a), then with 

i>(a) = (tpi(a)) and aW = (q^), (4.14) may be solved by the fixed point iteration scheme that 

can be described as 

«( fc ) = A -1 ^*" 1 )). (4.14) 


5 Numerical Examples 

In this section we present a numerical example using least-squares and interpolation to approx- 
imate k(s,t). Let k(s,t) = e st and f(t) = 1 - Then, the computed errors for the least 

squares method are shown in the following table. 



Errors 

n 

non-iterated 

iterated 

3 

.00622970884 

1.57623196E-05 

4 

.00356820367 

4.97616686E-06 

5 

.00230562003 

2.01950876E-06 

convergence rate sa 

2 

4 
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PROJECT 4 


HIDEAKI KANEKO, RICHARD D. NOREN and PETER A. PADILLA 

Superconvergence of Collocation Method for Hammerstein Equations 

The Hammerstein equation arises as a reformulation of a class of boundary value problems 
with nonlinear boundary conditions. The collocation method is one of the widely used numerical 
methods to approximate the solution of such equations due to its reasonable computation cost. 

We proved that the iterates of the collocation solutions for Hammerstein equation converge 
faster than the original collocation solution, a phenomenon commonly known as a supercon- 
vergence. This result extends the results obtained by Kaneko and Xu ( Superconvergence of the 
Iterated Galerkin Mathods for Hammerstein Equations -SIAM Jl. Num. Anal. June 1996 (to 
appear)) concerning a similar outcome for the Galerkin method. The degree of improvement in 
the rate of convergence depends upon the smoothness of the kernel involved. 
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SUPERCONVERGENCE OF THE ITERATED COLLOCATION METHODS 

FOR HAMMERSTEIN EQUATIONS 

H. Kaneko, R. D. Noren and P. A. Padilla 
Abstract 

In this paper, we analyse the iterated collocation method for Hammerstein equations with 
smooth and weakly singular kernels. The paper expands the study which began in [14] con- 
cerning the superconvergence of the iterated Galerkin method for Hammerstein equations. 

We obtain in this paper a similar superconvergence result for the iterated collocation method 
for Hammerstein equations. We also discuss the discrete collocation method for weakly sin- 
gular Hammerstein equations. Some discrete collocation methods for Hammerstein equations 
with smooth kernels were given previously in [3] and [18]. 

Key words: the iterated collocation method, the discrete collocation method, Hammerstein 
equations with weakly singular kernels, super convergence. 

Mathematics Subject Classification (1990): 65B05, 45L10. 

1 Introduction 

In this paper, we investigate the superconvergence property of the iterated collocation method for 
Hammerstein equations. In the recent paper [14], the superconvergence of the iterated Galerkin 
method for Hammerstein equations with smooth as well as weakly singular kernels was estab- 
lished. The paper generalizes the previously reported results on the superconvergence of the 
iterated Galerkin method for the Fredholm integral equations of the second kind [8], [9] [20]. A 
more important contribution made in [14] lies in the fact that the superconvergence result was 
established under weaker assumptions (Theorem 3.3 [14]). The approach used in [14] to establish 
the superconvergence of the iterated Galerkin method can easily be adopted to prove the results 
of Graham, Joe and Sloan [8], Joe [9] and Sloan [20] under weaker conditions imposed upon 
the Fredholm equations. This will be demonstrated in Section 3. In Section 2, we review the 
collocation method for Hammerstein equations as well as some necessary known results that will 
be pertinent to the matreials in the ensuing sections. We recall that the collocation method for 
weakly singular Hammerstein equations was discussed and some superconvergence results of the 
numerical solutions at the collocation points were discovered by Kaneko, Noren and Xu in [11]. In 
Section 3, the supererconvergence of the iterated collocation method for Hammerstein equations 
is established. The results obtained there encompass Hammerstein equations with smooth as well 
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as weakly singular kernels. Finally, in Section 4, we discuss the discrete collocation method for 
TT a.mm prst.ein equations with weakly singular kernels. The result obtained in this section extends 
the results of [3] and [18] which deals with the discrete collocation methods for Hammeratein 
equations with smooth kernels. Some examples are also included in this section. 

We note that there have been several other research papers published in recent years that 
describe various numerical methods for Hammerstein equations. A variant of Nystom method 
was proposed by Lardy [19]. The degenerate kernel method was studied by Kaneko and Xu 
■ [16]. We point out that a superconvergence of the iterates of the degenerate kernel method was 
recently observed when a decomposition of the kernel is done properly. This will be reported in 
a future paper [17]. The reader who is interested in more information on numerical methods for 
a wider class of nonlinear integral equations may find necessary materials in [2] and [5]. 

2 The Collocation Method 

In this section, the collocation method for Hammerstein equations is presented. Some materials 
from the approximation theory are also reviewed in this section to make the present paper self- 
contained. We consider the following Hammerstein equation 

x(t) - f k(t,s)ip(s,x($))ds = f(t), 0 < t < 1, (2.1) 

Jo 

where k , / and ip are known functions and x is the function to be determined. Define k t (s) = 
k(t , s) for t,s £ [0, 1] to be the t section of k. We assume throughout this paper unless stated 
otherwise, the following conditions on k, f and ip: 

1. lim t _* T || kt - MI°o = 0, r€ [0, 1]; 

2. M = sup t /o |fc(t,s)|ds < oo; 

3. / e C[0, 1]; 

4. ip(s, x) is continuous in s G [0, 1] and Lipschitz continuous in x G (— oo, oo), i.e., there exists 
a constant C\ > 0 for which 

|V>(s,xi) - ^(s,x 2 )| < Ci|xi - x 2 |, for all Xi,x 2 € (- 00 , 00 ); 

5. the partial derivative tp^ 0,1 ^ of ip with respect to the second variable exists and is Lipschitz 
continuous, i.e., there exists a constant Ci > 0 such that 

\ip( 0,1 \t,xi) — ip( 0,1 Xt,x 2 )\ < C 2 \xi — x 2 \, for all xi, x 2 € (- 00 , 00 ); (2.2) 
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6. for x E C[0, 1], ^(., x(.)), *(•)) € C[0, !]• 


We let 


(K9)(x)(t) = / k(t,s)i/j(s,x(s))ds. 

Jo 

With this notation, equation (2.1) takes the following operator form 


x - KVx = /. (2.3) 

For any positive integer n, we let 

n„ : 0 = to < *1 < ... < tn-l <tn = 1 

be a partition of [0,1]. Let r and v be nonnegative integers satisfying 0 < v < r. Let 5^(II n ) 
denote the space of splines of order r, continuity v, with knots at II n , that is 


Sr (n n ) = {x e C'fo, 1] : x|[ t . it . +1 ] € Vr- 1 , for each i = 0, 1, . . ., n - 1} 

where V T -\ denotes the space of polynomials of degree < r — 1. For the collocation method, we 
are interested in the cases v — 0 or 1. That is, it is posible to work with the space of piecewise 
polynomials with no continuity at the knots or with the space of continuous piecewise polynomials 
with no continuity requirement on the derivatives at the knots. We assume that the sequence of 
partitions II n of [0, 1] satisfies the condition that there exists a constant C > 0, independent of 
n, with the property: 

maxi <,-<„( t,- - L-i) < for ^ n ( 2 . 4 ) 

nun^, <„(*,• - L-i ) 

In many cases, equation (2.1) possesses multiple solutions (see e.g. [16]). Hence, it is assumed 
for the remainder of this paper that we treat an isolated solution xo of (2.1). Let /,• = (t,_i,t,) 
for each i = 1, . . . , n. Then for v = 0, we let r,i,r, 2 , . . .,Ti r be the Gaussian points (the zeros of 
the rth degree Legendre polynomial on [—1,1]) shifted to the interval Ii. We define 

<^0 = {rif 1 < * < n, 1 < j < r}. (2.5) 

The points in Go give rise to the piecewise collocation method where no continuity between 
polynomials is assumed. This is the approach taken by Graham, Joe and Sloan [8]. Joe [9], on 
the other hand, considered the continuous piecewise polynomial collocation method. His method 
corresponds with taking v = 1. Here we define the set G\ of the collocation points to be the 
set consisting of the knots along with the Labatto points (the zeros of the first derivative of the 
r — 1th degree Legendre polynomial) shifted to the interval Namely, let £ r _i = 1 and for 
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1 < Z < r — 2 (r > 3), let £/ denotes the Zth Labatto point. If |J,j denotes the length of then 
G\ contains 

T (i-i)(r-i)+i+i = ^(*.-1 + U + 1 < * < n, 1 < / < t - 1, and r x = t 0 = 0 (2.6) 

The analyses of [8] and [9] are very s imil ar. We therefore confine ourselves to developing the 
collocation method for Hammerstein equations that is analogous to the method of [8]. An 
obvious extension to the continuous piecewise collocation method will be left to the reader. 
Define the interpolatory projection P n from C[0, 1] + 5^ (II n ) to 5^ (II n ) by requiring that, for 

x € c[o, i] + s*(n„), 

P n x(Tij) = x(Tij), for all Tij £ Go. (2.7) 

Then we have, for x £ C[0, 1] + 5^(II n ) 

P n x -* x, as n — ► oo (2.8a) 

and consequently 

sup||P n || < c. (2.86) 

n 

The collocation equation corresponding to (2.3) can be written as 

x n - P n Vx n = P n f (2.9) 

where x n £ S” (II n ). Now we let 

. Tx = / H- K^lx 

and 

T n x n = P n f + P n KVx n 

so that equations (2.3) and (2.9) can be written respectively as x = Tx and x n = T n x n . We 
obtain; 

Theorem 2.1 Let xq £ C[0, 1] be an isolated solution of equation (2.3). Assume that 1 is not an 
eigenvalue of the linear operator (K^)'(xq), where (K'9 ) , (xq) denotes the Frechet derivative of 
Kty at xo. Then the collocation approximation equation (2.9) has a unique solution x n £ B(xo,S) 
for some 6 > 0 and for sufficiently large n. Moreover, there exists a constant 0 < q < 1, 
independent of n, such that 

^ < ||*. - zolloo < (2-10) 

where a n = ||(7 - T4(x 0 )) -1 (T n (x 0 ) - T(xo))||oo- Finally, 

E n (x 0 ) < ||x n - xolloo < CE n (x 0 ), (2.11) 

where C is a constant independent of n and E n (x 0 ) = inf ue x„ |jxo — ti||oo- 
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A proof is a straight application of Theorem 2 of Vainikko [23] and is demonstrated in the 
proof of Theorem 2.1 [11]. We denote by W™[ 0, 1], 1 < p < oo, the Sobolev space of functions g 
whose TO-th generalized derivative belongs to Z p [0, 1]. The space W™[ 0, 1] is equipped with 
the norm 

m 

llslkf = L ll^ll,- 

k=0 

It is known from Demko [6] and De Vore [7] that if 0 < f < r, 1 < p < oo, m > 0 and x G W™, 
then for each n> 1, there exists u n G 5^(II n ) such that 

Ik-«n||p<c^||x||^, (2.12) 

where p = min{m,r} and h = maxi <,<„(£,■ — The inequality (2.12) when combined with 
Theorem 2.1 yields the following theorem; 

Theorem 2.2 Let xq be an isolated solution of equation (2.3) and let x n be the solution of 
equation (2.9) in a neighborhood of x o. Assume that 1 is not an eigenvalue of (K^l)\x 0 ). If 
x 0 € W, Iq, then 

11*0 ~ *n||oo =.0(/i M ), 

where \i — min{/, r}. If aj 0 € Wj, (1 < p < oo), then 


||*0 ~ *n||oo = 0{h v ), 

where v — min{/ — 1, r}. 

When the kernel k is of weakly singular type, namely if 

k(t,s) = m(t,s)g a (\t- s\), (2.13) 


where m G C M+1 ([0, 1] x [0, 1]) and 

9 a {s ) 


a®" 1 , 0 < a < 1, 
logs, a = 1. 


(2.14) 


then the solution xo of equation (2.3) does not, in general, belong to W™. To better characterize 
the regularity of the soution of (2.3) with weakly singular kernel, consider a finite set S in [0, 1] 
and define the function u>s(t) = inf{|t — s| : s G S}. A function x is said to be of Type(a,fc,5), 
for — 1 < a < 0, if 

|* W (*)I < CM*)] 0 "* t i 5, 

and for a > 0, if the above condition holds and x G Lip(a). Here Lip(a) = {x: \x(t) - *(s)| < 
C|t — s| a }. It was proved by Kaneko, Noren and Xu [12] that if / is of Type((3,p,{ 0,1}), then 
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a solution of equation (1.1) is of Type( 7 ,//, {0,1}), where 7 = min{a, (3}. The optimal rate of 
convergence of the collocation solution x n to xq can be recovered by selecting the knots that are 
defined by 

U = (l/ 2 )( 2 i/n) ? , 0 < i < n/ 2 , (2J5) 

ti-l - t n -i, n /2 <i<n , 

where q = r / 7 denotes the index of singularity. Details can be found in [11]. 

3 The Iterated Collocation Method 

The faster convergence of the iterated Galerkin method for the Fredholm integral equations of 
the second kind compared to the Galerkin method was first observed by Sloan in [21] and [22]. 
On the other hand, the superconvergence of the iterated collocation method was studied in [ 8 ] 
and [9]. Given the equation of the second kind 

x — Kx = /, (3.1) 


where K is a compact operator on X = C[ 0, 1] and x, / € X, the collocation approximation x n 
is the solution of the following projection equation 

x n - P n Kx n = P n f. (3.2) 

Here P n is the interpolatory projection of (2.7). The iterated collocation method obtains a 
solution x' n by 

< = f + Kx n . (3.3) 


Under the assumption of 

\\KP n — K\\ — ► 0 as n — *• oc (3.4) 

it can be shown that 

||x - x'J < ||(/ - KP n )-'\\\\K{x - P„x)||. (3.5) 

The assumption (3.4) is satisfied if X = L 2 and P n is the orthogonal projection satisfying 
|| P n g — < 7 1 1 — *■ 0 for all g in the closure of the range of the adjoint K * of K since in this case 
||A'P n - K\\ = || P n K* - K*\\. The results of Sloan were recently generalized to the iterated 
Galerkin method for Hammerstein equations by Kaneko and Xu [14]. The main theorem of [14], 
Theorem 3.3, that guarantees the superconvergence of the iterates was proved by making use of 
the collectively compact operator theory. 
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The purpose of this section is to study the superconvergence of the iterated collocation 
method. For the collocation solution x n of (2.9), we define 

x’ n = f + K*x n . (3.6) 

A standard argument shows that x' n satisfies 

x' n = f + KVP n x' n . (3.7 a) 

We denote the right side of (3.7a) by S n x' n , namely 

S n x' n = f + K^P n x' n . (3.7 b) 

We recall the following two lemmas from [14]. 

Lemma 3.1 Let xo G C[ 0, 1] be an isolated solution of (2.3). Assume that 1 is not an eigenvalue 
of (K^)'(xq). Then for sufficiently large n, the operators I — S' n (xo) are invertible and there exists 
a constant L > 0 such that 

\\(I -S'^xoT'Woo < L, for sufficeintly large n. 

Lemma 3.2 Let x 0 G C[0, 1] be an isolated solution of equation (2.3) and x n be the unique 
solution of (2.9) in the sphere B(x 0 ,6i). Assume that 1 is not an eigenvalue of (K^)'(x 0 ). Then 
for sufficiently large n, x' n defined by the iterated scheme (3.6) is the unique solution of (3.7) in 
the sphere B(x 0 , S). Moreover, there exists a constant 0 < q < 1, independent of n, such that 

lTg - - Xollo ° - r^? 

where (3 n - ||(/ - 5^(a:o)) _1 [5„(xo) - r(® 0 )]||oo- Finally, 

IK - zolloo < CE n (x 0 ). 


The definitions of 6 and are described in [14]. Following the development made in [14], we 


ij>(s, y) = -ip(s, y 0 ) + ^ (0ll) (s, y 0 + 9{y - y 0 ))(y - y 0 ), 


where 0 := 0(s,yo,y ) with 0 < 9 < 1. Also let 


9(t, s, 2/o, V , 0) = k(t, s)^ (0,1) («, 2/o + 0(y ~ 2/o)), 

(G n x)(t)= f g(t,s,P n x 0 (s),P n x' n (s),9)x(s)ds, 

Jo 

and ( Gx)(t ) = J 0 a g t (s)x(s)ds, where g t (s) = k(t 1 s)ip^°' 1 \s,xo(s)). Now we are ready to state 
and prove our main theorem of this paper. The proof is a combination of the idea used in [14] 
(Theorem 3.3) and the one used in [8] (Theorem 4.2). 
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Theorem 3.3 Let x 0 € C[0, 1] be an isolated solution of equation (2.3) and x n be the unique 
solution of (2.9) in the sphere B(x 0 ,6i). Let x' n be defined by the iterated scheme (3.7). Assume 
that 1 is not an eigenvalue of (K9 )'(xq). Assume that xq € W[ (0 < / < 2 r) and g t € W™ 
(0 < m < r) with ||5t||vr 1 m bounded independently oft. Then 

||x 0 - x^lloo = G(/i 7 ), where 7 = min{l,r + m}. 

Proof: From equations (2.3) and (3.7), we obtain 

*0 - *n = K (* x o - ^n<) = K(*X 0 - VP n x 0 ) + K(VP n x 0 - *P n x' n ). (3.9) 

Using (3.8), the last term of (3.9) can be written as 

K(VP n x o - *P n x' n )(t) = (G n P n (x o - x' n ))(t). 

Equation (3.9) then becomes 

x 0 - x' n = K(Vx o - $P„x 0 ) + G n P n (x 0 - x' n ). (3.10) 

Using the Lipschitz condition (2.2) imposed on ^ 0,1 \ for x e C[0, 1], 

||(G n x) — (Gx)||oo < C 2 SUp f |fe(^,5)|rf5|jx|| 00 (||P n Xo — Io|too + H-fnllooll^n ~ ^o||oo)- 

0<t<l Jo 

This shows that 

|| Gn - Glloo < AfC 2 (||P„xo - Xolloo + c||x'„ - Xolloo) -*• 0 as n —y oo. 

Also, for each x € C[0, 1], 

sup |(GP„x)(t) — (Gx)(t)| = sup | f g t {s)[P n x{s) - x(s)]ds\ < MMi\\P n x - xWoo, 

0<t<l 0<t<l Jo 

where 

Mi = sup |^ o,1 )(t,x o (0)l < +oo. 

0 <t<l 

It follows that GP n — *■ G pointwise in C[0, 1] as n — *■ oo. Again since P n is uniformly bounded, 
we have for each x € C[0, 1], 

||G n P„X - Gxllco < ||G n - G||oo||Pn||oo||*||oo + ||GP n X - Gx|U. 

Thus, G n Pn — *• G pointwise in G[0, 1] as n — ► oo. By Assumptions 2, 5, and 6, we see that there 
exists a constant C > 0 such that for all n 

M 0ll) (s, PnXo(s) + 0{PnxM ~ P n X 0 {s ))) | < C 2 ||P n X 0 - X 0 ||oo + 0C 2 P\\x' n - X 0 ||oo + M X < C. 
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This im plies that {G„P n } is a family of collectively compact operators. Since G = (K^)'(x 0 ) is 
compact and (I - G)~ l exists, it follows from the theory of collectively compact operators that 
(/ - G n P n )~ l exists and is uniformly bounded for sufficiently large n. Now using (3.10), we see 
that 

lko-<||oc<C'||iir(*xo-fP„xo)li. 

Hence we need to estimate lli^xo - ®Pn®o)||- The following four inequalities are known (The- 
orem 4.2 [8]). Let r/} n € 5°(n n ) be such that 

E ||(*0 - tn^WwrW Z ch l -’\\x 0 \\ wl , o <j<l, (3.11) 

t=l 

mg* II V4 j) llws(/f) ^ c ll*°llw}, 3 ^ °- ( 3 - 12 ) 

Also for each t € [0, 1], there exists <p n<t € 5^(II n ) such that 


E Ufa* “ < c h m J K m , 0 <j<m, 

i = 1 

llrf t !llwg(/<) < cKm, 3 > 0, 

where K m = sup 0 < t<1 \\k t \\ w ™ < oo. Now for t € [0, 1] we have 


AT($X 0 - $P„Xo)(t) = ( gt - <Pn,t,X 0 - PnXo) + (V>n,t, (I ~ -Pn)(*0 - V>n)) 

+(<Pn,t,(I - PnWn)- 


(3.13) 

(3.14) 


(3.15) 


Using equations (3.11)-(3.14) along with the arguments from [8] (p.362) we can show that each 
of the three terms is bounded by ch? uniformly in t. This completes our proof. □ 

One way to establish the superconvergence of the iterated collocation method for the Fredholm 
equation is to assume (3.4). In the context of the present discussion, (3.4) is equivalent to 
assuming 

||(Jf$y(x 0 )(/ - Pn)|c[a,6]l|oo — 0 as 71 ♦ OO. (3.16) 


Theorem 3.3 was thus proved under weaker assumptions. The idea used to prove Theorem 3.3 
originates from [4] (section 6) in which the superconvergence of the iterated collocation method 
for the Fredholm equations was established by showing that {KP n } is a family of collectively 
compact operators. 

Finally in this section, we investigate the superconvergence of the iterated collocation method 
for weakly singular Hammerstein equation. Specifically, we consider equation (2.3) with kernel 
given by (2.13) and (2.14). An enhancement in the rate of convergence is given in the following 
theorem. 
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Theorem 3.4 Let xo £ C[ 0, 1] be an isolated solution of equation (2.3) and x n be the unique 
solution of (2.9) in the sphere B(xo, 6\) with kernel defined by (2.13) and (2.14) and knots defined 
by (2.15). Let x' n be defined by the iterated scheme (3.7). Assume that 1 is not an eigenvalue of 
(KV)'{x 0 ) and that ^ 0,1 )(-,2:o(-)) is of Type(a,r, {0,1}) for a > 0 whenever xo is of the same 
type. Then 

Proof: We follow the proof of Theorem 3.3 exactly the same way to (3.15), which is 

K(Vx 0 - VP n X 0 )(t) = ( g t - <p n ,t, Xo - P n x o) + (<Pn,t, ( I ~ P n )( x 0 - tin)) 

+(<Pn,U (I ~ PnW n). 

The difference in superconvergence arises from the degree to which we may bound the first term. 
As in Kaneko and Xu [14] (Theorem 3.6), using an argument similar to [15], it can be proved 
that there exists u € S^(U n ) with knots II„ given by (2.15) such that \\g t — u|ji = 0{h a ). Here 
h = maxi< t < n {* t - — Xi-i}. Then 

\(gt ~ Tnji X Q ~ Pn x o)\ ^ ||ff< ~ ^n.tlllll-^O — Pn^olloo 

= 0(h a+r ). 

The rest of proof follows in the same way as described in [8] (p.362).D 

4 The Discrete Collocation Method for Weakly Singular Ham- 
merstein Equations 

Several papers have been written on the subject of the discrete collocation method. Joe [10] 
gave an analysis of discrete collocation method for second kind Fredholm integral equations. A 
discrete collocation- type method for Hammerstein equations was described by Kumar in [18]. 
Most recently Atkinson and Flores [3] put together the general analysis of the discrete colloca- 
tion methods for nonlinear integral equations. In this section, we describe a discrete collocation 
method for weakly singular Hammerstein equations. In the aforementioned papers [10, 18, 3], 
their discussions are primarily concerned with integral equations with smooth kernels. Even 
though, in principle; an analysis for the discrete collocation method for weakly singular Ham- 
merstein equations is similar to the one given in [3], we feel that a detailed discussion on some 
specific points pertinent to weakly singular equations, -e.g.,a selection of a particular quadrature 
scheme and a convergence analysis etc, will be of great interest to practioners. Our convergence 
analysis of the discrete collocation method presented in this section is different from the one given 
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in [3] in that it is based upon theorem 2 of Vainikko [23]. The idea of the quadrature used here 
was recently developed by Kaneko and Xu [15] and a complete Fortran program based on the 
idea is being developed by Kaneko and Padilla [13]. A particular case of the quadrature schemes 
developed in [14] is concerned with an approximation of the integral 

1(f) = f f(x)dx, (4.1) 

Jo 

where / € Type(a,2r,S ) with a > —1. For simplicity of demonstration, we assume 5 = {0}. 
We define q = and a partition 

it a :x Q = 0,xi = n~ q ,xj = j 9 h, j = 2,3,...,ra. (4.2) 


Now we construct a piecewise polynomial S T of degree r — 1 by the following rule; S r (x) = 0, 
x £ [x 0 ,ii) and S r (x) is the Lagrange polynomial of degree r - 1 interpolating / at {u^}j =1 for 
x € [®i,x, + i), i = 1, 2, . . ., n — 2 and for x € [x„_i,x n ]- Here denote the zeros of the 

rth degree Legendre polynomial transformed into [x,-,x t+1 ). Our approximation process consists 
of two stages. First, 1(f) is approximated by 

1(f) = / f(x)dx = f(x)dx. (4.3) 

Jxi i=1 Jxi 

Second, 1(f) is approximated by I(S T ) = S r (x)dx. A computation of I(S T ) can be accom- 
plished as follows; let s: [x,-,x,- + i] —»■ [—1, 1] be defined by s = 2x so that 

/(/) = /_//(»)* (4.4) 


where 

n_1 1 1 1 
F f( s ) = 53 2^ i+1 “ + 2^ i+1 + 

t=l 

If {it,-: i = 1, 2 . . . , r) denotes the zeros of the Legendre polynomial of degree r, then 

«=1 

with U(x) the fundamental Lagrange polynomial of degree r — 1 so that 

- -A fi , 

I(S r ) = YwiFf(ui), where = I lfx)dx. (4.5) 

i=i *' -1 

It was proved in [15] that 

W/)-l(5r)l = 0(™- 2 '). (4.6) 
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In this section, we examine equation (2.1) with the kernel k defined by (2.13) and (2.14). When 
the knots are selected according to (2.15), as stated earlier, it was shown in [11] that the solution 
x n of the collocation equation (2.9) converges to the solution x of (2.1) in the rate that is optimal 
to the degree of polynomials used. Specifically, x n must be found by solving 

* n («W) “ / SaflUj 0 ~ s\)m(uf\s)ip(s,x n (s))ds = (4.7) 

JO 

where i = 0, 1, .. .n — 1 and j = 1,2, .. .r. 

The discrete collocation method for equation (2.1) is obtained when the integral in (4.7) is 
replaced by a numerical quadrature given in (4.5). Let fc,j(s) = g a (\u^ — s|)m(u^,s). Then 

fo S , «(l u i ) ~ s\)™( u j >, s)4>(s, x n (s))ds = kij(s)ip(s, x n (s))ds 

, (4.8) 

= fo J + /*< 0 x n {s))ds. 

u i 

The integrals in the last expression of (4.8) represent two weakly singular integrals which can be 
approximated to within 0(n~ 2r ) order of accuracy by (4.5) by transforming them to [—1, 1] and 
selecting the points in (4.2) appropriately. 

Writing (4.7) as 

P n Xn - PnKVx n = P n f, (4.9) 

we consider the approximation x n to x n defined as the solution of 

®n = QnX n — PnK n ^lX n + Pnf \ (4.10) 

where K n is the discrete collocation approximation to the integrals in (4.8) described above. 

We will use Theorem 2 of [23] to find a unique solution to (4.10) in some 6 neighborhood of x n , 
where n is sufficiently large. Clearly, Q' n (x) = P n K n ^'(x), where $'(a;)[y](s) = ip 0,l \s, x(s))y(s). 
For sufficiently large n, (4.9) has a unique solution in some 6 neighborhood of x. To see that 
I - Q' n (x n ) is continuously invertible with {(/ — ^(zn)) -1 }^# uniformly bounded, it is enough 
to observe that {Q^^n)}^! is collectively compact, and to do this we will show that 

I QnM[x](t) - Q' n (x n )[x){t') | = | P n K n 9'(Xn)x(t) - P n K n *' (Xn)z(t') j- 0 (4.11) 

as t — *■ f, for each x € C[ 0, 1], [1]. Here N is some sufficiently large number. 

If we show (4.11), then part (a) of Theorem 2 [23] is also verified. In order to verify part 
(b) of Theorem 2 [23], we only need to establish (because of the uniform boundedness of {(/ — 

<5n(a: n )) _1 }^N) that 

II Qn( x ) — Qni. x n ) llooS: L || X — X n ||oo$ LS, (4.12) 
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for some constant X, and 

II Qn(x n ) - T n {x n ) ||-» 0 as n — oo. (4.13) 

Once this is done, Theorem 2 [23] applies yielding a unique solution x„ in some neighborhood 
of x n (for sufficiently large n) and 

|| x n - x n ||< La n < X || Q n (x n ) - T n (x n ). H* . (4.14) 

(Here and throughout the remainder of the section, X denotes a generic constant, the exact 
value of which may differ at each occurance.) This inequality will be used to obtain the order of 
convergence. 

Considering (4.11), the right hand side is bounded by T\ + + T 3 , where 

T\ =| P n K n 9'(x n )x(t) - P n KV'(x n )x(t) I, 

T 2 ^P n K*Xx n )x{t)-P n K*\x n )x(!?)\, - 

T 3 =| P n K n V'(x n )x(t') - P n K n *'(x n )x{U) | . 

Let e > 0. Since {P n }£Li is uniformly bounded, T\ + T 3 < ^ by applying (4.6) with f(s) = 
x n (s))x(s) and letting n be sufficiently large. For T 2 we have 
T 2 <M /J I k(t,s) - ktfjs) | ds < M(Si + S 2 ), 
where 

Si = f ft*(| * - 1 1) | m(t, s ) - m(t', s) | ds 
Jo 

and 

S 2 = / | Scd t - s |) - g a { 1 1' - s |) || m(t',s ) | ds. 

Jo 

but 

Si < sup | m(t, s) — m(t', s) \ I g a (\t — s |)ds 
o<s<i Jo 

< X sup | m(t,s ) — m(t',s ) |— *■ 0 as t — ► t', 

0 < 5<1 

and 

S 2 <LJo \g a (\t-s\)-g a (\lf -s\)\ds 

= £{i t a - (t'r 1 + 1 a - tr - (1 - ?r \+£\t-t n 

— ► 0 as t — * t' . 

Hence (4.11) holds. For (4.12), 

|| Q' n (x) - Q'n(x n ) ||oo = || PnK n (*'(*) ~ *'(*n)) ||< MC \\ X - X n ||< M6 = q < 1 
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for 6 s uffi ciently small. Note that we have used the uniform boundedness of {P n },{hf„} and 
because ^ (°’ 1 )(s, y(s)) is locally Lipschitz, so is the operator 

W : C[0, 1] — ► B(C[ 0, 1],C[0, 1]) (the space of bounded linear operators from C[0, 1] into C[0, 1]). 
For (4.13), we have 

|| Qn(x n )-T n (x n ) ||oo=|| P n (K n 9x n - K9x n ) ||< L II (K n -K)9(x n ) < L(R X + R 2 + R 3 ) (4.15) 
where 

R x =|| K n *(x n )-K n *(x 0 ) ||, R 2 =|| K n V(x 0 )-K9(x 0 ) ||, Rz =|| KV(x 0 )-KV(x n ) || . (4.16) 
But 

Rx<L || *(x«) - *(x 0 ) ||< CiL || x n - so II (4.17) 

because W is a Lipschitz operator and {K n } is uniformly bounded, and also 

R 3 < M II «(x 0 ) - «(*n) ||< C X M II x n - II - (4.18) 

Finally, 

R 2 = 0(n" 2r ) (4.19) 

by (4.6) using f(s) = $(x,x 0 (s)). 

Thus Vainniikko’s Theorem yields a unique solution x n for n sufficiently large and (4.14) 
holds. Now (4.14) and (4.15) - (4.19) show that 

||x n -x n ||=0(R-^) (4.20) 


where (3 is the minimum of 2r and the order of convergence of || xo — x n ||. We summerize the 
results obtained above in the following theorem: 

Theorem 4.1 Let xo be an isolated solution of equation (2.3) and let x n be the solution of 
equation (2.9) in a neighborhood of x o- Moreover, let x n be the solution of (4-10). Assume that 
1 is not an eigenvalue of (K^)'(x o). If xq € then 


11*0 - *n||oo = 0(h M ), 

where p = min{/,r}. If x 0 (E Wj, (1 < p < oo), then 

||*0 ~ *n||oo = 0{h v ), 


where v — min{/ — 1, r}. 
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PROJECT 5 


HIDEAKI KANEKO, RICHARD D. NOREN and PETER A. PADILLA 

Singularity Preserving Galerkin Method for Hammerstein Equations With Logarithmic Ker- 


in the recent paper Singularity preserving Galerkin methods for weakly singular Fredholm 
integral equations, Jl. Int. Eqs. and Appl. 6 (1994) SOS-334, Y. Cao and Y. Xu established the 
Galerkin method for weakly singular Fredholm integral equations that preserves the singularity 
of the solution. Their Galerkin method provides a nu m erical solution that is a linear combination 
of a certain class of basis functions which includes elements that reflect the singularity of the 
solution. The purpose of this paper is to extend the result of Cao and Xu and to establish 
singularity preserving Galerkin method for Hammerstein equations with logarithmic kernel. 
First, a singularity expansion for the soultion of Hammerstein equation with logarit hmi c kernel 
is given. Secondly, this singularity expansion is used to obtain the numerical Galerkin scheme 
that preserves the singularity of the solution. 

An application is given to a Dirichlet problem with a certain class of nonlinear boundary 
conditions. Numerical experiments are being performed. 


1 



Singularity Preserving Galerkin Method For 
Hammerstein Equations With Logarithmic Kernel 

Hideaki Kaneko* Richard D. Noren 
Department of Mathematics and Statistics 
Old Dominion University 
Norfolk, Virginia 23529-0077, USA 
and 

Peter A. Padilla 

NASA Langley Research Center 
MS-130 

Hampton, Virginia 23681-0001, USA 


‘This author is partially supported by NASA under grant NCCl-213 


0 



SINGULARITY PRESERVING GALERKIN METHOD FOR HAMMERSTEIN 
EQUATIONS WITH LOGARITHMIC KERNEL 


H. Kaneko, R. D. Noren and P. A. Padilla 
Abstract 

In a recent paper [3], Y. Cao and Y. Xu established the Galerkin method for weakly singu- 
lar Fredholm integral equations that preserves the singularity of the solution. Their Galerkin 
method provides a numerical solution that is a linear combination of a certain class of basis 
functions which includes elements that reflect the singularity of the solution. The purpose 
of this paper is to extend the result of Cao and Xu and to establish singularity preserving 
Galerkin method for Hammerstein equations with logarithmic kernel. An application is given 
to a Dirichlet problem with nonlinear boundary condition. 
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1 Introduction 


In this paper, we are concerned with the problem of obtaining a numerical solution of weakly 
singular Hammerstein equations with logarithmic kernel by the Galerkin method that preserves 
the singularity of the exact solution. Namely we establish a method that generates an approxi- 
mate solution in terms of a collection of basis functions some of which are comprised of singular 
elements that reflect the characteristics of the singularity of the exact solution. The idea of the 
method originates in the recent paper by Cao and Xu [3]. Cao and Xu studied the characteristics 
of the singularities that are pertinent to solutions of the weakly singular Fredholm equations of 
the second kind. Let C[0, 1] denote the space of all continuous functions defined on [0, 1]. The 
weakly singular Fredholm integral equations of the second kind can be described as 

y(s)~ f 9 a{\s - t\)m(s,t)y(t)dt = f(s), 0 < s < 1 (1.1) 

Jo 

where / £ C[0, 1], m is sufficiently smooth and 


g a (\s-t\) = 


|s - f|“, -1 < a < 0, 

log|s-<|, a = 0, 


( 1 - 2 ) 


where y is of course the function to be determined. It is well documented (see, e.g. [16], [13], [4], [20]) 
that the solutions of the equations described in (1.1) exhibit, in general, mild singularities even 
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in the case of a smooth forcing term /. Here by “mild” singularities, we mean the singularities 
in derivatives. The papers of Richter [13] and Graham [4] contain singularity expansions of the 
solutions of equation (1.1) in the case of m(s,t) = 1. The results of Graham were recently 
generalized by Cao and Xu for equation (1.1). Information concerning the type of singularities 
that solutions have is useful when solving equation (1.1) numerically. In order to approximate 
functions with mild singularities, many investigators utilized the important theorem of Rice [14] 
that gives an optimal order of approximation to such functions. Rice’s theorem is that of non- 
linear approximation by splines using variable knots. Based upon this idea of approximating 
the solutions of equation (1.1) by splines defined on nonuniform knots, the collocation method, 
the Galerkin method and the product- integration method were established for equation (1.1) by 
Vainikko and Uba [21], by Graham [4] and by Schneider [17] respectively. A modified collocation 
method was introduced in [12] which also uses the idea of Rice. Recently there has been some 
considerable interest in the study of the following weakly singular Hammerstein equation: 

y(s)~ [ g a (\s-t\)m(s,t)ip(t,y(t))dt = /(s), 0 < s < 1 (1.3) 

Jo 

where /, m and g a are defined as in (1.1) and (1.2) and ip is a known function. We will see in 
Section 4 that equation similar to (1.3) arises naturally in connection with Dirichlet problem with 
certain nonlinear boundary conditions. A study on the regularities of the solution y of equation 
(1.3) is reported in [9], extending the results of [16]. Subsequently, Kaneko, Noren and Xu 
used the regularity results to establish the collocation method for weakly singular Hammerstein 
equations in [10]. The approximate solutions provided by these methods for equations (1.1) 
and (1.3) are in the form of piecewise polynomials that are not always satisfactory as a tool 
for approximating functions with singularities. This observation is quite evident in the areas of 
finite element analysis. Hughes and Akin [6] list several problems (e.g. ‘upwind’ finite elements 
for treating convection operators [5], [8], [7]; boundary-layer elements [1] etc.) in which the finite 
element shape functions are constructed to include polynomials as well as singular functions. 
Singular shape functions are introduced to the set of basis functions through asymptotic analysis 
on the solution of the problem that is being considered. It should be pointed out that the 
analysis involved in the aforementioned papers on the finite element method is centered around 
the collocation method. The problems such as the choice for the extra collocation points for 
singular basis elements or the rate of convergence are not addressed in these papers. It should 
be pointed out that the location of additional collocation points for singular basis elements is 
critical in detemining the rate of convergence of numerical solutions. A detailed discussion on 
this subject can be found in [11]. A singularity preserving collocation method, because of the 
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reasons mentioned above, seems to be more difficult to establish. 

In this paper, a singularity expansion for the solution of equation (1.3) with logarithmic 
kernel is given. This extends the results in [9] and [3]. Only the logarithmic kernel is considered 
here because of our primary interest in obtaining numerical solution of a Dirichlet problem 
with nonlinear boundary condition as described in Section 4. It is a routine matter, however, 
to establish, following the argument of Section 2, a singularity expansion for the solution of. 
equation (1.3) with algebraic, singularity. A detail is left to the reader. The paper is organized, 
as follows: in Section 2, we study the regularity property of the soultion of (1.3) and establish 
its singularity expansion. The results obtained there generalize the results of [3] and [9]. The 
singularity expansion is then utilized in Section 3 to achieve the singularity preserving Galerkin 
method for equation (1.3). Finally, in Section 4, the singularity preserving Galerkin method is 
applied to a class of Dirichlet problems with nonlinear boundary condition. Examples are also 
included in this section. 

2 Singularity Expansion for Weakly Singular Hammerstein Equa- 
tions 

In this section, we consider the following Hammerstein equation with logarithmic singularity, 

y(s) - f log |s - t\m(s, y(t))dt = f(s), 0 < s < 1 (2.1) 

Jo 

(see (1.3) also). We let 

K^y(s) = f log | s - f |m(s, y(t))dt. (2.2) 

Jo 

Then equation (2.1) can be written in operator form as 

y - KVy = /. (2.3) 

Let H n denote the Sobolev space H n [ 0, 1] = {w : € L2[0, 1]} equipped with the norm 

IMIff" - (nui ^ where describes the zth generalized derivative of w. We also 

let W = W n be the linear space spanned by the functions s' log- 7 s, (1 - s)’iog- 7 (l - s);i,j = 
1,2, ...,n— 1. Throughout the paper, we assume the following conditions: 

m e C 2n ([0, 1] x [0, 1]), n > 1, m e C a ([0, 1] x [0, 1]), n = 0. (2.4) 

€ C 2n+1 (R x R) (2.5) 
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few®H n . 


( 2 . 6 ) 


We define 

Ky(s)= f log - t\m(s,t)y(t)dt. (2.7) 

Jo 

Also let tij(s) = s p log 9 s, and ti 2 (s) = (1 - s) p log 9 (l — s), where p, q > 1 axe integers. First we 
quote the following result (lemma 4.4(2)) from [3]. 

Lemma 2.1 Let / £ H n ~ x and assume m £ C n+1 ([0, 1] X [0, 1]). Then, 

(Kf)( s ) = Y [c? 5,7 lo S * + dj( 1 - s) 3 log(l - s)] + u„(s), 
j = i 


n— 1 9+1 n— 1 

(A'ui)(s)= Y ]Ec,jS J (logs)' + Y d ii l ~ S Y lo g(l “ s ) + «#(*)» 

j=p+i i=i i=9+i 


and 


n— 1 


n—l q+l 

(ku 2 )( s )= e E c ii( 1 - s ) i ( lo s(i-*)) i + E djS 3 log s -1- u n (s). 

j= P +i i=i j=,+i 


Lemma 2.2 If ui(s) — s p log g s, U 2 (s) = (1 - $) r log u (l — s), where p,q,r,u > 1 are integers, 
then uiu 2 £ W © H n . 

Proof: Expand Ui in series about t = 1 and u 2 about t = 0 : 

«i (o = Er^o 1 6 *( 1 - + a (o» “2co = e"=o «•■<*’ + hit), 

= Pi(t) + /i(0 — P 2 (0 + h{t) 

where f[ k \t) = 0((1 — t) n ~ k ) near t = 1,/j is analytic at t = 1, and ~ tt^OO - P^( 0) as 

t — *■ 0+; f 2 k \t) = 0(t n-A: ) near t = 0,f 2 is analytic at t = 0, and f 2 k \t) ~ tt^CO — -P 2 ^(0 35 
t-l - . 

Now tii ti 2 = P 1 P 2 + P 1/2 + P 2 /i + /i/ 2 . Clearly PiP 2 is in £T n . For /i/ 2 , we have 
Each term /i^(0/ 2 ” *^(0> 1 0 “ satisfies 

/i ( %)/| n ~ 0 (O = 0(A (,) (OP) = O([4 i} (0 - A (0 (0)]p) - 0 

a s<-+ 0+. 

Similarly 
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/i' ) (*)/2 n "' ) (0 - o as t -> 1-. Thus f x f 2 £ C n C P n . For / a P 2 we have fi(t)P 2 (t) = 
(ui(t) — P x (t))P 2 (t) = Ui(t)P 2 (t) — Pi{t)P 2 (t). Since P 2 is a polynomial, u x € W, it is easy to 
see that uiP 2 € W © H n (see [[3], (4.7)]). So f\P 2 £ H n . Similarly f 2 P\ £ W © //”, and Lemma 
2 has been verified. □ 


Lemma 2.3 A product of an H n function with a function in W is in H n © W . 


Proof: Let g £ H n and let u\ and u 2 be defined as before prior to Lemma 1. For gu x we write 

ui(s)g(s) = E^o 1 am°l s i+r log * s + /• g H ((r )( s _ a )n-i da 

= T X + T 2 . 


Since 7i £ W © P n , we turn to T 2 and write 

= HU © &r\* log" <][(» - 1)).. Jt] £ S <">(<7)(. - 
+s p log 9 sg( n )(s). 

But s p log ? s £ L°°,g^ £ X 2 [0, 1] so (s p log 9 s)g^ n \s) £ L 2 . 

For the terms 

K(s) = ^[s p log* s]jf g^ n \o)(s - of^do 
we have, for some constant M and nonegative integer a 

I Mo) I 

= Ms(- logs)”* /„■ | I dcr. 

But g( n ) £ L 2 [0, 1], so by Hardy’s inequality [15] (p. 72) j/q | g^(er) | do £ L 2 [0, 1]. Since 
s(— logs) 0 £ L°° it follows that b n £ i 2 [0, 1]. Hence £ L 2 [0, 1], or T 2 £ H n , This proves 
gu x £ W © H n . 

The case for gu 2 £ W © H n is similar. O 

Finally we need the following: 


Lemma 2.4 The operator K^l maps W © H n into W © H n+1 . 


Proof: Let y = w + h, w £ W, h £ H n . We use Taylor’s theorem in the form 

H x ) = S i^ (fc) (®)(* ~ a) k + -7 / (* - tr ) n ^ n+ 1 )(<r)d<7. ( 2 . 8 ) 

fc=o K - n ' Ja 


5 



Letting x = y(s) and a = h(s ) allows us to write • 


(jr«)(»)(0 = U=o n fo 1 log I < - * I m(t, s)^ k )(h(s))w(s) k ds 

+^i fo log I t - s I m(t, s) n+1) (a)(y(s ) - a) n dods (2.9) 

= ELo£ Mt) + ±B(t). 

By (3), i/>W(h(s)) € H n ,k = 0, 1, n, and by expanding w(s) k with the multinomial expan- 
sion, it is clear that w(s) k is a sum of terms in W as well as terms of the form as p log 9 5(1— . 
s) r log u (l — s), p,q,r,u > 1 are integers. The constant, a, depends on p,q,r and u. Since 
^l fc l(h(s)) € H n and w(s) k 6 W ® H n ,k = 0,1, it follows from Lemma 3 that 

^ k \h(s))w(s) k eW®H n . (2.10) 


By Lemma 1 and (2.10), we have 


A k eW®H n+1 . 


For B(t), if we prove that 


( 2 . 11 ) 


F(s) = r (s) ^ n+1 \(T)(y(s) - o) n do G W © H n , 

Jk(s) 

then, also by Lemma 1, B(t) = I([F](t) will be in IF © H n+1 . This will complete the proof of 
this lemma. First of all, suppose n > 1. We write 


F'(s) = -^ n+1 \h(s))w{s) n h'(s). 

Since h G H n ,ip G C 2n+1 , rp( n+1 \h(s)) G H n . By Lemmas 2 and 3, — ip( n+1 \h(s))w(s) n G 
H n ® W. Since h! G H n ~ x , it follows that — ^ n+1 \h(s))w(s) n h'(s) G H n ~ l ®W (Lemma 
2), Since F' G H n ~ l © W it is clear that F G H n © W. Second, let n = 0. Then F(s) = 
/ a(/j ^'(cr)da = tp{y{s)) - ip{h(s)) € L 2 [0, 1 ] C W © H°. 

Thus 

B(t) eW®H n . (2.12) 

By (2.9), (2.11) and (2.12), it follows that K^l maps W © E n into W © U"* 1 . □ 

Using the lemmas which we proved above, we obtain the following main result of this section. 

Theorem 2.5 Suppose the conditions (2.4)-(2.6) hold and y is an isolated solution of (2.1). 
Then there are constants aij and b{j, for i,j — 1,2, ...,n— 1, and there is a function v n in H n 
such that 

n—1 n— 1 

y(t) - Y lo 6 J t + bij( 1 - ty log J (l - t )] + u n (0- (2.13) 

;=i j = i 
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Proof: For n = 0, this follows from Lemma 4 with n — 0. Assume that the result holds 

for n = k, that is, if / G H k © W*, then (2.13) holds with n = k. Say y = Wk + Vk, where 

v k e = Ei=i & l Kf W * + M 1 - W(1 - 01- 

Now consider the case n = /; + 1 and suppose / € ® W'fc+i- 

Since y = Wk + Vk we write y = 5^2/ + / = K^!(wk + «*) + /• From Lemma 1, /f + £ 

VFfc+i ® H k+1 . The proof is complete. □ 

3 Singularity Preserving Galerkin Method 


In this section, we establish the singularity, preserving Galerkin method for equation (2.1). First 
we recall the definition of the space of spline functions of order n. Define the partition of [0, 1] as 

A : 0 = to < ti < ... < tk < tfc+x = 1. 


Let 


h = max (ti — L_i), 
i<;<fc+i v ' 

and assume h —* 0 as k — ► oo. Denote by II n the set of polynomials of degree n — 1. Then the 
space of splines of order » with knots ti ' s of multiplicity n — 1 — v is defined as 

Si = sp* (A) = {s e C"[0, 1] : s| fi € n n , }, 


where 0 < v < n — 1 and /,• = {ti-\,t%) for i = l,2,...,fc+ 1. It is well known that the dimension 
of S% is d = n{k + 1) — k( 1 + 1 /). is spanned by a basis consisting of 5-splines {5j}f =1 . We let 


V£ = W®SZ (3.1) 

and denote the orthogonal projection of £ 2 ( 0 , 1] into V£ by P^. The singularity preserving 
Galerkin method for approximating the solution of equation (2.3) requires the solution yh € V£ 
that satisfies the following equation: 

y h - P£KVy h = iff. (3.2) 

More specifically, we need to find yh in the form 

n— 1 n— 1 d 

yh(s) = Y lo S J s + J2 /M 1 - • s )‘ lo g J (! - *) + Y 7«' 5 «( s ) ( 3 - 3 ) 

*.i=i »,i=l »=i 
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where {a tJ , Pij}*jli an <l are found by solving the following system of nonlinear equations: 

E” 7=1 log J s p log’ a) + E- Jil Ai((l - sy log j (l - s), s p log’ s) 

1 <*«;*' fog J 5 + E"jii fti(! - 5 )‘ fog^ 1 - s)), s p log? s) 

= (/,s p log’s) P,9= l,2,...,n- 1 

E"jil Q t'i("S* log J 5, (1 - s) p log’(l - s)) + Erjii Aj((l - «)*' fog J (l - 5), (1 - s) p log’(l - 5)) 
-(^^(E”jii a ii 5 ‘ fog J 5 + E"jii Pi A 1 - 5)*’fog J (l - «)),(! - s) p log’(l - a)) 

= (/,(l-s) p log’(l -«)) p,q = 1,2 1 

E"7ii <*iAj fog J *, B p ) + E-Ji, Ai((l - 5)*' log j (l - s), fi p ) 

-(**(E?jii fog J * + E,” 7=i M* - *) f fog J (i - «)), 5 P ) 

= (f,B p ) p = l,2,...,d 

where (•, •) denotes the usual inner product defined on L 2 [ 0, 1]. Now let Ph be the orthogonal 
projection of L 2 [0, 1] into S£. Then we have 

PhV — ► v for all v £ L 2 [ 0, 1]. (3.4) 

It is well known (e.g. [18]) that if g € H n , n > 0, then for each h> 0, there exists fa € S% such 
that 

||5 ~ <f>hh < Ch n \\g\\H n , (3.5) 

where C > 0 is a constant independent of h. By virtue of the fact that P^u is the best L 2 
approximation of u from S£, we see immediately that 

|| PhV ~ tt||2 < ||u - <^|| 2 < Ch n ||uj| H n, for all u € H n . (3.6) 

The following lemma from [3] is useful in the sequel. 

Lemma 3.1 Let X be a Banach space. Suppose that U\ and U 2 are two subspaces of X with 
U\ C U 2 . Assume that Pi : X —+ U\ and P 2 : X — *■ U 2 are linear operators. If P 2 is a projection, 
then 

II* “ -P2*|U < (1 + ||F 2 ||x)||* - A*|U for all x ex. 

For convenience, we introduce operators T and T n by letting 

Ty = f + KVy 

and 

T h y h = if/ + P?KVy h 


(3.7) 

(3.8) 
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so that equations (2.1) and (3.2) can be written respectively as y = Ty and yh = T^yh- The 
following theorem guarantees the existence of a solution of the singularity preserving Galerkin 
method (3.2) and describes the accuracy of its approximation. 

Theorem 3.2 Let y € L 2 [0, 1] be an isolated solution of equation (2.1). Assume that 1 is not 
an eigenvalue of the linear operator (K9Y(y), where (K^)'(y) denotes the Frechet derivative of 
K $ at y. Then the singularity preserving Galerkin approximation equation (3.2) has a unique 
solution yh, such that ||j/ — yh\\ 2 < 6 for some 6 > 0 and for all 0 < h < ho for some ho > 0. 
Moreover, there exists a constant 0 < q < 1, independent of h, such that 

S-'Sh-nhig? <»> 

where ah = ||(/-T^(j/)) -1 (T 1 f l (?/) - T(y)) \\ 2 . Finally, ify = w + v with w € W and v € H n , then 

||y — 2^112 < C/i n ||t;||iyn, whenever 0 < h < ho, (3.10) 

where C > 0 is a constant independent of h. 

Proof: The existence of a unique solution yh of equation (3.2) and the inequalities in (3.7) can 
be proved using Theorem 2 of Vainikko [19]. A detailed argument can be found in [10]. To get 
(3.10), first we note from Lemma 3.1, for v € L 2 [0, 1], 

Now, from (3.9), 

\\y - yfcih 

= T^\\(I-n(y))-\Th(y)-f(y))\U 
< C\\PfK*y - KVy + iff - f\\ 2 
= G\\Pfy - y\\ 2 . 

where C is independent of h. Using the uniform boundedness of {P®}, (3.11) and (3.12), we 
obtain 

Ib-Jftlla < Ch n \\v\\ Hn . 

□ 


(3.11) 


(3.12) 


4 Nonlinear Boundary Value Problem 


In this section, we consider the following Laplace’s equation with nonlinear boundary condition 
in R 2 : 

A u(P) =0, Pf O 

£t(P) =-WP,u(P)) + f(P), 
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Pen 
p e r = dD 


(4.1) 



where D is a simply connected open region in R 2 with an open contour T and np denotes the 
exterior unit normal to T at P. The functions / and ij) are given and we assume that / 6 C(r) 
and ij) € C(r x R). The solution u is to be found in C 2 (D) H C X (D). The function tj) is assumed 
to be continuous on T x R. The problem (4.1) in the case of a closed smooth boundary T was 
considered by Atkinson and Chandler in [2]. They employed the method of piecewise polynomial 
product integration and that of trigonometric product integration to approximate the solution 
of (4.1). In the current problem in which T describes a boundary that is open, one expects 
logarithmic singularities in u at the two ends of the boundary. Now it is well known that, using. 
Green’s representation formula for harmonic functions, the function u satisfies 

“< p > = h l - 5 . i |p - «*«» (4 - 2) 

for all P € D. Using the boundary condition in (4.1) and letting P approach to a point of T, we 
obtain 

<‘IQ)j^ioe\P-QlMQ)-iJri’l,QMQ))log\P-Q\<b{Q) 

= -iJrf(Q) l °&\P-Q\MQ), Per. 

We denote the double layer operator by T, 

Tv(P)=:^J^v(Q)-~log\P-Q\da(Q), Pel, (4.4) 

and the single layer operator by 5, 

Sv(P) = -~ / v{Q)log\P-Q\d<j(Q), Per. (4.5) 

If we put $v(P) = if>(P,v(P )), P € T, then equation (4.3) can be written in operator form as 

u - Tu + SV{u) = S /. (4.6) 


Define the parametrization r(t) = (f(t), Tj(t)) for t € [0, 1] and assume that r € C°°[0, 1] and 
|r'(t)| / 0 for t € [0, 1]. The double layer operator and the single layer operator now become 
respectively, 


1 f 1 

Tv(t) = i / 
7T Jo 




[£(«) - f(01 2 + Ms) - *?(01 2 

for v € C[0, 1] with the kernel having the value 


v(s)ds 


(4.7) 


{'(*)<(*) -l/(t)g"(Q 

mw+vm 


when s = t and 

Sv(t) = — — f log |s — t|v(s)|r , (s)|ds. 
7T Jo 


(4.8) 


10 



The kernel of T is well behaved whereas the kernel of S needs our scrutiny. As in [2], we write 
Sv as 


Sv(t) = - £ fg v(s)|r'(s)|{log)t — <s| + log 1 1 — s + f| + log 1 1 — t + s|}ds 

Jo U ( 5 M S )I 1°S [ |,_ t |( l 1 l^[) ( (x ) Lt+,) ] ds - 

Due to the fact that the double layer operator T contains infinitely differential kernel, using the 
arguments to prove Lemma 2.4, we see that T + S9 maps W ® H n into W ® H n+1 . Hence we 
obtain the following result that parallels the results obtained in Theorem 2.5.. 


Theorem 4.1 Suppose the conditions (2.5) and (2.6) hold and u is an isolated solution of (4-6). 
Then there are constants aij and b(j, for i,j = 1,2, ..., n — 1, and there is a function v n in H n 
such that 

n— 1 n— 1 


y(t) = 2 Efci* 1 lo g J * + M 1 - O' i°g J (i - 0] + Mt)- C 4 -^) 

*=i j= i 


Define 

t u = Tu - SV(u) + Sf (4.11) 


so that equation (4.6) can be written as 


u — Tu. 


(4.12) 


The singularity preserving Galerkin method is now described by 

u h - PffTu h + P?S*(u h ) = P h Sf. (4.13) 

where Uh € Vf) and P® is the orthogonal projection of I/2[0, 1] into Vjf as defined in the previous 
section. By letting 

f h u = P?Tu-lfSV(u) + P?Sf, (4.14) 

equation (4.13) can be written as 

Uh = ThUh- (4-15) 

The following theorem guarantees the existence of the solution of equation (4.15) and describes 
its accuracy as an approximation to u that is the soution of equation (4.12). A proof is an easy 
exercise of modifying the argument given in the proof of theorem 2.1 [10], hence we omit it. 

Theorem 4.2 Let u € L2[0, 1] be an isolated solution of equation (4-6). Assume that 1 is not 
an eigenvalue of the linear operator T + (S^)'(u). Then the singularity preserving Galerkin 
approximation equation (4-13) has a unique solution uh such that ||u — Uh \\2 < 6 for some 6 > 0 
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and for all 0 < h < h 0 for some ho > 0. Moreover, there exists a constant 0 < q < 1, independent 
of h, such that 

< Il«-«ft||2 < — (4.16) 
where ah = ||(J — f' h (u))~ 1 (Th{y) — T(y))|| 2 . Finally, if y = w + v with w £ W and v € H n , then 

||tt - tifclh < Ch n \\v\\nn t whenever 0 < h < h 0 , (4.17) 

where C > 0 is a constant independent of h. 
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